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SUMMARY 


A  theoretical  and  experimental  study  of  correlation  bandwidth  over 
troposcatter  paths  was  made  over  a  seventeen  month  period  from  February 
1969  through  June  1970.  The  field  test  program  consisted  of  experiments 
to  determine  cross  correlation  coefficient  versus  frequency  separation 
and  fade  characteristics  over  four  different  troposcatter  paths  in  the 
northwestern  region  of  New  York  State.  X-and  C-band  signals  were  trans¬ 
mitted  simultaneously  using  a  common  antenna  from  the  RADC  Test  Site  in 
Model  City,  New  York,  to  Ontario  Center,  New  York,  to  Whitford  Field  near 
Weedsport,  New  York,  to  Point  Petre,  Ontario  Canada  (over  water),  or  to 
Port  Byron,  New  York.  Each  of  the  paths  were  of  different  characterisi tcs 
so  that  the  effect  of  terrain  and  path  length  could  be  assessed.  The 
receiving  instrumentation  was  located  in  a  mobile  trailer  and  moved  from 
site  to  site  every  several  weeks  to  provide  seasonal  variations.  The 
resulting  data  were  recorded  on  magnetic  tape  for  computer  processing  at 
the  Martin  Marietta  computer  facility  in  Orlando,  Florida.  The  output 
information  consisted  of  simultaneous  X-and  C-band  measurements  of  cross 
correlation  coefficient  versus  frequency  separations  up  to  9  MHz  and 
statistics  of  fade  rate,  fade  duration,  fade  depth,  and  signal  amplitude 
distributions . 

The  data  from  the  field  test  program  showed  that  the  correlation  band¬ 
width  of  the  X-band  signals  and  the  C-band  signals  were  about  the  same. 
Therefore,  the  same  order  of  frequency  diversity  is  obtainable  in  the 
same  bandwidth  at  the  higher  frequency  if  10  foot  antennas  are  used  over 
these  paths  which  are  typical  of  the  MALLARD  troposcatter  links.  There 
is,  of  course,  the  expected  additional  path  loss  and  the  medium- to- aper¬ 
ture  coupling  loss  of  the  associated  antennas.  The  fade  rates,  durations, 
etc.  ,  were  also  studied.  It  was  found  that  the  fade  rates  were  sometimes 
enormous  and  beyond  the  measurement  equipment  limitation  of  about  60  Hz. 

As  expected,  the  X-band  fade  rates  exceeded  the  C-band  rates  in  propor¬ 
tion  to  their  frequencies.  The  durations  and  other  statistics  are  best 
described  as  typical  of  what  would  be  expected  from  existing  theory. 

In  the  theoretical  study  a  literature  search  was  made  and  a  theoretical 
model  for  the  prediction  of  correlation  bandwidth  in  terms  of  identifiable 
parameters  was  constructed.  Several  models  were  studied.  It  was  determined 
that  all  models  appear  to  agree  reasonably  well  if  a  suitable  scale  factor 
is  included.  A  Gaussian  model  was  constructed  of  the  form: 

P  =  exp[-  (Sf  A  )  2] 

where:  p  =  the  envelope  cross  correlation  coefficient  for  a 

frequency  displacement  of  f  Hertz 

A  =  the  multipath  spread  in  seconds 

S  = 


and 


a  constant. 


r 


This  particular  model  is  convenient  because  it  uses  the  simple  exponential 
form  that  can  be  handled  mathematically  in  other  equations  where  p  is  a 
factor.  The  term  A  is  calculated  from  a  computer  program  that  traces  the 
upper  and  lower  rays  of  the  troposcatter  effective  beamwidth  through  the 
atmosphere  to  determine  the  difference  in  the  two  paths  and,  hence,  the 
multipath  spread.  The  effective  beamwidth  was  assumed  to  be  up  to  one  half 
the  scatter  angle  in  the  original  model. 

A  comparison  of  data  from  the  field  test  program  showed  that  S  was  not 
originally  a  constant,  but  varied  considerably  from  one  path  to  another. 
Further  study  indicated  that  the  effective  beamwidth  used  in  the  ray  tracing 
computer  program  was  itself  a  function  of  the  scatter  angle  rather  than 
simply  one  half  of  it.  Using  the  field  test  data  obtained  from  the  four 
paths,  an  empirical  function  was  determined  from  which  the  effective  beam- 
width  could  be  calculated.  This  function  when  applied  to  the  ray  tracing 
technique  to  determine  multipath  spread  resulted  in  a  constant  value  for 
S  in  the  model.  This  empirical  model  yields  good  agreement  with  the 
observed  data. 

It  was  particularly  noted  that  the  correlation  coefficient  versus  fre¬ 
quency  separation  was  stationary  only  for  short  periods.  From  one  15  minute 
test  to  the  next  the  measured  correlation  bandwidth  may  move  through  factors 
of  easily  2  to  1.  From  this  it  is  postulated  that  two  mechanisms  cause 
scattering;  one  of  these  is  turbulent  and  the  other  is  a  more  stable  lay¬ 
ering  phenomenon.  Furthermore,  it  is  believed  that  varying  combinations 
of  the  two  scattering  mechanisms  are  actually  the  result  of  what  is  seen 
at  the  receiver.  The  narrowest  correlation  bandwidths  appear  to  be  due  to 
turbulent  scattering  in  an  unstable  common  volume,  while  the  wider  corre¬ 
lation  bandwidths  might  be  attributed  primarily  to  the  layering  phenomenon. 

For  comparison,  some  of  the  correlation  coefficient  data  obtained 
from  a  previous  MALLARD  troposcatter  program  is  included  herein.  This 
added  data  is  from  the  ECOM  link  between  the  Hexagon  and  Tobyhanna, 
Pennsylvania,  which  is  a  low  scatter  angle ,  100  mile  path. 

The  appendices  contain  the  Rochester  weather  data  obtained  for  the 
period  of  testing.  Also  included  is  the  correlation  bandwidth  computer 
program  listing  which  is  written  in  FORTRAN  for  the  IBM  360  computer.  A 
study  of  the  Bello  Model  for  correlation  bandwidth  is  included  and  applied 
to  these  paths.  It  was  found  that  the  Bello  Model  did  not  fit  the  observed 
data  very  closely  although  it  appears  to  be  a  very  good  theoretical  develop¬ 
ment.  A  complete  test  log  is  included  for  reference  purposes. 
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I. 


PROGRAM  OBJECTIVES 


Project  MALLARD  is  a  program  for  the  design  and  development  of  a 
digital,  automatically  switched  communication  system  for  the  field  Army 
with  characteristics  as  described  in  the  project  MALLARD  Quadripartite, 
Communications  Plan,  and  Proposals  for  Research  and  Development.  This 
digital  system  will  be  designed  to  transmit  full  duplex  voice,  teletype, 
facsimile,  and  digital  data  at  various  speeds  with  total  traffic  security. 
The  communication  links  between  points  will  use  many  different  devices 
such  as  cable,  microwave,  radio,  etc.,  as  appropriate  for  the  particular 
link.  Troposcatter  paths  up  to  200  km  nominal  are  to  be  used  in  areas  where 
direct  cable  or  simple  microwave  systems  are  not  applicable. 

The  objective  of  this  program  was  to  obtain  data  necessary  for  the 
design  of  troposcatter  modems  for  use  in  future  MALLARD  troposcatter 
systems.  In  support  of  this  objective,  propagation  data  were  collected 
over  four  types  of  troposcatter  paths  to  empirically  determine  cross 
correlation  versus  frequency  spacing  and  other  fade  statistics.  These 
data  provide  a  basis  for  determining  the  maximum  bit  rates  that  can  be: 
satisfactorily  transmitted  by  frequency  diversity  methods  with  a  stated 
error  probability  under  known  conditions  of  path  length,  path  terrain, 
season,  antenna  size  and  beamwidth,  radio  frequency  band,  occupied  band¬ 
width,  radiated  power,  and  frequency  spacing.  These  data  are  also  directly 
applicable  to  the  prediction  of  space  diversity  maximum  bit  rates  at 
stated  error  probabilities  since  the  correlation  bandwidth  and  fade  sta¬ 
tistics  directly  affect  the  error  probabilities  and  bit  rates  obtainable. 

The  propagation  data  were  collected  through  the  use  of  a  fixed  C- 
and  X-band  transmitting  system  located  at  the  RADC  troposcatter  test 
site  at  Model  City,  New  York.  The  receiving  instrumentation  was  located 
in  a  special  van  that  was  moved  to  various  locations  to  provide  variations 
in  terrain,  path  length,  and  season. 
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II  ANALYTICAL  MODEL  FOR  THE  CORRELATION  BANDWIDTH 
A.  SUMMARY  AND  CONCLUSIONS 

As  part  of  this  program,  an  empirical  model  for  the  correlation  band¬ 
width*  has  been  developed  based  on  the  subsequently  presented  experimental 
data.  As  an  aside,  this  development  has  entailed  a  study  and  comparison 
with  generally  accepted  models  of  the  troposcatter  channel ,  particularly 
those  proposed  by  Booker-Gordon^ ,  Friis-Crawford-Hogg2 ,  Rice2,  Sunde4,  and 
Bello5 . 

The  model  proposed  here,  being  empirical,  is  believed  valid  for  paths 
similar  to  those  included  in  the  field  tests,  that  is,  paths  in  the  range 
140  to  200  km  (85  -  125  statute  miles).  For  longer  paths  an  extrapolation 
is  possible  and  the  model's  validity  has  been  ascertained  for  a  few  cases 
for  which  data  were  available. 

One  problem  in  predicting  correlation  bandwidth  that  has  been  insuffi¬ 
ciently  addressed  in  the  past  lies  in  the  marked  time  variation  in  the  cor¬ 
relation  function  for  a  specific  link.  The  correlation  bandwidth  is  inversely 
related  to  the  multipath  spread  and  experimental  evidence  points  to  the  fact 
that  this  spread  is  normally  not  constant  for  more  than  a  few  hours.  Over 
a  fixed  path  this  variation  must  be  due  to  changes  in  meteorological  condi¬ 
tions,  particularly  in  the  common  volume.  At  times  it  appears  that  the 
scatterers  in  the  common  volume  are  in  turbulent  motion  so  that  the  entire 
geometrical  common  volume  produces  effective  scattering.  At  other  times 
there  is  evidence  of  strongly  refracting  horizontal  layers  restricting  the 
size  of  the  common  volume  and  consequently  the  multipath  spread.  Turbulent 
scattering  would  then  result  in  a  correlation  bandwidth  small  in  comparison 
with  the  case  of  stratified  layers. 

It  has  been  experimentally  determined  that  the  correlation  function  can 
be  modeled  as  Gaussian  with  variance  differing  by  a  factor  of  about  2.5  for 
the  two  extremes  of  turbulent  and  layer  scattering.  Thus  the  approach  has 
been  taken  of  obtaining  the  correlation  function  for  the  turbulent  case  with 
an  adjustment  of  the  variance  for.  the  stratified  case.  The  median  correla¬ 
tion  bandwidth  falls  between  these  two  bounds.  It  is  argued  that  turbulent 
scattering  corresponds  experimentally  to  the  observed  lower  deciles  while 
layering  predominated  in  the  upper  deciles. 

When  no  pronounced  layers  are  present,  the  behavior  of  the  link  can  be 
accounted  for  by  changes  in  the  gradient  of  the  radio  refractive  index. 

♦Defined  here  as  the  frequency  separation  at  which  the  envelope  correlation 

p  is  0.4. 
e 
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This  point  has  been  investigated  in  detail  and  computed  tables  that  allow 
estimation  of  the  frequency  correlation  function  for  a  proposed  link  on  the 
basis  of  the  surface  refractive  index  Ns  have  been  investigated.  Further¬ 
more,  equations  are  given  that  allow  Ns  to  be  computed  from  readily  avail¬ 
able  meteorological  data  such  as  temperature,  pressure,  and  relative  humid¬ 
ity.  The  presence  of  layers  is  not  easily  predicted  from  surface  values  of 
Ns  and  will  depend  on  sonde  data. 

In  addition  to  affecting  the  take-off  angles  at  the  receiver  and  trans¬ 
mitter  sights,  the  terrain  is  important  in  its  effect  on  the  gradient  of 
radio  refractive  index.  Since  this  gradient  can  be  estimated  from  Ns,  the 
terrain  becomes  important  in  its  bearing  on  local  weather  conditions  and 
therefore  Ng . 

The  scatter  angle  has  been  found  to  be  the  single  most  important  param¬ 
eter  in  predicting  correlation  bandwidth  lower  deciles.  The  actual  antenna 
beamwidth  is  believed  to  have  only  a  minor  significance  and  in  its  place  the 
concept  of  effective  beamwidth  has  been  introduced.  It  is  defined  in  terms 
of  the  scatter  angle.  For  practical  links,  the  effective  beamwidth  is  gen¬ 
erally  smaller  than  the  actual  half  power  beamwidth. 

B.  CORRELATION  BANDWIDTH  EXPRESSIONS  OBTAINED  FROM 
PREVIOUS  SCATTERING  MODELS 

The  operation  of  troposcatter  radio  communication  systems  depends  on 
electromagnetic  energy  scattered  by  the  earth's  troposphere  due  to  inhomo¬ 
geneities  in  its  dielectric  constant.  The  frequencies  used  for  transmission 
are  typically  in  the  1  to  8  GHz  range.  For  these  frequencies,  it  is  rela¬ 
tively  easy  to  construct  narrow  beam  antennas.  Figure  1  shows  a  typical 
troposcatter  link  with  all  relevant  parameters  defined.  In  Figure  1  it  has 
been  assumed  that  the  antennas  are  located  at  ground  level  and  that  there 
are  no  obstructions  along  the  horizon.  The  center  of  the  antenna  beams 
would  normally  be  slightly  elevated  above  the  horizon. 

Active  theoretical  investigation  of  the  troposcatter  mechanism  started 
in  the  late  1940 's  with  work  by  Booker  and  Gordon'*'.  They  postulated  a  model 
involving  isotropic  turbulence  in  the  common  volume  as  the  scattering  medium. 
Based  on  this  model  they  and  others  developed  theories  of  expected  channel 
behavior . 


Figure  1.  Troposcatter  Link  Geometry 
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A  convincing  model  was  subsequently  developed  by  Friis,  Crawford  and 
Hogg  -  in  which  horizontal  layers  with  sharp  transitions  in  the  gradient  of 
the  dielectric  constant  were  postulated.  This  model  is  realistic  in  that 
the  •i'-rity  of  wind  induced  air  currents  will  tend  to  flow  parallel  to  the 
e?ct  •  '  s  rface  and  thus  often  lead  to  horizontal  layer  formation. 

ic  is  to  be  expected  that  the  Booker  and  Gordon  model  will  apply  when 
the  troposphere  is  well  mixed,  as  it  may  be  on  a  stammer  afternoon,  while 
the  Friis,  Crawford  and  Hogg  model  will  apply  when  the  troposphere  is 
strongly  stratified.  This  condition  is  more  likely  to  occur  during  the 
winter.  The  first  interim  report®  gave  a  comparison  of  the  Booker-Gordon 
and  Friis-Crawford-Hogg  scattering  theories . 

The  correlation  bandwidth  of  a  fading  radio  channel  is  a  measure  of 
the  usable  bandwidth  of  the  channel.  It  also  indicates  how  far  apart  in-band 
frequency  diversity  channels  must  be  spaced  to  obtain  substantial  diversity 
gain. 


The  troposcatter  channel  suffers  from  two  kinds  of  fluctuations.  One 
fluctuation  is  rapid  fading  caused  by  the  random  motion  of  the  scatterers. 

This  fading  has  a  typical  time  scale  in  fractions  of  seconds.  Another  fluc¬ 
tuation  is  a  slow  variation  that  follows  variations  in  weather,  day,  and 
season.  The  averaging  time  for  measurement  of  the  frequency  correlation 
function  must  be  chosen  to  include  many  rapid  fades,  but  not  so  great  as  to 
cover  a  substantial  period  of  the  slow  fluctuations.  A  range  of  suitable 
averaging  times  is  from  2  to  10  minutes. 

The  frequency  correlation  function  is  related  mathematically  to  a  more 
easily  visualized  function,  the  multipath  spread.  The  multipath  spread  is 
the  time  elongation  an  ideal  impulse  would  suffer  if  transmitted  over  the 
channel.  If  the  multipath  spread  for  the  channel  can  be  determined,  the 
correlation  bandwidth  can  be  calculated  from  it.  As  shown  in  Figure  1,  it 
appears  to  be  a  simple  matter  to  determine  the  maximum  multipath  spread. 

It  is  only  necessary  to  evaluate  the  difference  in  propagation  delay  between 
signals  arriving  via  the  top  and  bottom  of  the  common  volume.  However,  this 
approach  proves  quite  inaccurate  since  radio  waves  do  not  travel  in  straight 
lines  in  the  troposphere.  Furthermore,  since  the  power  transfer  is  a  sensi¬ 
tive  function  of  the  scatter  angle,  the  effective  common  volume  will  in  gen¬ 
eral  be  smaller  than  the  geometrical  one.  It  appears  that  the  multipath  spread 
is  a  function  of  both  the  path  length  and  scattering  angle,  but  evidence 
gained  here  and  elsewhere  indicates  that  it  does  not  depend  very  strongly  on 
antenna  beamwidth.  Experimental  results  have  shown  only  a  small  change  in 
correlation  bandwidth  of  a  4  GHz  link  as  the  receiving  antenna  diameter  was 
changed  from  8  to  CO  feet^. 

In  the  theoretical  models  developed  by  S.  0.  Rice^  and  also  by  Hirai, 
Fukushima,  and  Kurihara®  it  was  assumed  that  the  scatterers  were  randomly 
distributed  throughout  the  common  volume  according  to  a  three  dimensional 
Gaussian  probability  law.  Therefore,  this  model  takes  into  account  the 
effects  of  antenna  beamwidth  but  does  not  take  into  account  the  dependence 
of  scattered  power  on  scattering  angle.  Nor  does  it  take  into  account  the 
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effect  of  horizontal  layers  in  the  common  volume.  Nevertheless,  with  an 
appropriate  definition  of  the  common  volume,  the  model  does  agree  with 
other  theoretical  models  as  well  as  with  the  results  of  experiments. 

According  to  their  models,  the  envelope  correlation  coefficient  between 
two  frequencies  separated  by  f  in  Hz  is : 


p  (f) 
e 


-  (2  Ttcrf ) 2 
e 


(i) 


The  critical  parameter  in  equation  (1)  is  a,  which  depends  on  the 
maximum  difference  in  propagation  delay  from  the  various  parts  of  the 
effective  common  volume.  The  effective  common  volume  has  been  shown  ex¬ 
perimentally  to  depend  mainly  on  path  length  and  scattering  angle.  But 
because  of  refraction  of  the  antenna  beams  by  the  earth's  atmosphere, 
simple  straight  line  geometrical  reasoning  will  usually  not  produce  very 
accurate  results. 


A  model  has  been  developed  for  the  common  volume  that  will  take  this 
refraction  into  account.  A  simplified  version  of  this  model  uses  an  effec¬ 
tive  earth  radius 
this  report. 

_  2 1  .  e 

C  Sln  2 


R  ,  which  is  defined  and  derived  in  another  section  of 
0 

In  Rice's  model  the  multipath  spread  is  given  by 


a 


(2) 


where 


t~  radius  of  the  common  volume/  VT 

0  =  =  scattering  angle 

Re 

c  =  velocity  of  light  in  vacuum 


d  =  distance  between  transmitter  and  receiver. 


Since  equation 


(2)  can  be  simplified  as 


cr 


(3) 


From  Figure  1  it  can  be  seen  that  the  radius  of  the  effective  common  volume 
is  approximately  ad/4 ,  where  a  is  the  effective  beamwidth  of  the  antennas. 
Thus , 


and 
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Because  of  the  scattering  mechanism  the  effective  antenna  beamwidth 
a  often  is  only  weakly  related  to  the  free  space  beamwidth  of  the  antennas 
in  use.  For  multipath  purposes,  a  workable  expression  for  the  effective 
antenna  beamwidth  given  by  Gordon  9  is: 

“  -  I  •  <6> 


This  relation  is  only  an  approximation  and  will  be  modified  later  in  this 
section.  However,  when  this  value  is  substituted  for  a  in  equation  (5)  we 
obtain: 


d  9 
8  c  VT  * 


The  Rice  Gaussian  model  will  now  be  compared  with  one  developed  by 
Sunde4.  Sunde  also  assumed  independent  scattering  from  a  large  number  of 
scatterers  within  the  common  volume.  By  use  of  a  Fourier  series  type  of 
argument ,  Sunde  obtained : 

pe,f)  -  (8) 

where  A  is  the  maximum  departure  from  the  average  transmission  delay. 


A  general  expression  can  be  derived  for  the  maximum  differential  trans 
mission  delay,  based  on  geometrical  considerations.  This  equation  is  valid 
even  when  the  transmitting  and  receiving  antennas  have  different  beamwidths 
The  notation  used  is  shown  in  Figure  2. 


If 


8  =  total  differential  path  length 


smin  =  minimum  path  length 
smax  =  maximum  path  length 


Figure  2.  Troposcatter  Path  With  Unequal  Antennas 
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then  it  is  easily  shown  (Ref.  6)  that 

2 

smin  =  a  Z  d(l  +  |  )  ,  9^1,  (9) 

cos  9  8 

2 

and 

smax  *  d(l  -  1/2  (|  +  aR)  (|+aT)),  (10) 


so  that  the  differential  path  length  is 

8  *  |smax  -  smin|-  f  [ <f  +  °R>  t|  +  “t>  '  |  ]•  <U> 

Equation  11  gives  the  maximum  difference  in  path  length  even  when  the 
antennas  used  at  each  end  of  the  path  have  different  beamwidths. 


The  time  dispersion  relative  to  the  average  time  delay  for  the  case  of 
transmitting  and  receiving  antennas  of  unequal  beamwidths  is  then  approxi¬ 
mately  given  by. 
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If  it  is  somewhat  arbitrarily  assumed  that  cr=  0.573A,  then  Rice's 
expression  for  the  correlation  function  becomes 


P  (f)  =  exp  [  -  (2  7T  '  0.573  fA)2] 
e 


(14) 


and  under  this  assumption,  the  Rice  and  Sunde  (Equation  8)  models  are  in 
excellent  agreement  as  shown  in  Figure  3.  The  two  curves  are  almost  iden¬ 
tical  over  the  region  of  most  interest,  pe  >0.4.  To  see  how  well  the  two 
theories  agree,  the  values  obtained  for  A  and  o’  must  be  compared.  From 
equations  7  and  13,  it  is  noted  that: 


A 

rr 


3  d  92 
16  c 

d  A  2 


=  \/T~  =  2.58 


while  under  the  assumption  that  cr  =  0.57  3 ,  the  result  obtained  is: 


A  =  1 

<7  0.573 


1.75  . 


Thus,  although  the  two  mathematical  models  were  obtained  by  quite  different 
methods,  it  is  significant  that  they  are  in  reasonably  good^uantitative 
agreement.  A  more  detailed  comparison  of  the  S.  O.  Rice  and  E.  D.  Sunde 
developments  for  the  envelope  frequency  correlation  function  was  presented 
in  the  first  interim  report6. 

A  third  method  for  predicting  the  frequency  correlation  function  for 
a  troposcatter  link  is  to  assume  that  simple  geometric  ray  tracing  is  cor¬ 
rect  and  actually  integrate  over  the  common  volume  formed  by  the  intersec¬ 
tion  of  the  transmitting  and  receiving  antenna  beams.  This  has  been  done 
among  others  by  Bello6.  These  results  are  compared  (for  a  200  mile  tropo¬ 
scatter  link  operating  at  a  carrier  frequency  of  850  KHz  and  using  28  foot 
antennas)  with  the  predictions  of  the  Rice  model  in  Figure  4.  The  general 
shape  of  the  curve  obtained  by  integrating  the  scattered  power  over  the 
common  volume  is  very  similar  to  the  one  obtained  by  Rice,  but  it  is  appar¬ 
ent  that  the  scale  factor  is  different.  The  Bello  model  is  discussed  in 
more  detail  in  Appendix  A  as  it  pertains  to  the  paths  used  in  this  program. 


Figure  3.  Comparison  Between  Rice's 
and  Sunde' s  Model  for  the  Envelope 
Frequency  Correlation  Function 


Figure  4.  Comparison  Between  Bello's 
and  Rice’s  Model  for  the  Envelope 
Frequency  Correlation  Function 


In  conclusion,  it  is  to  be  noted  that  all  three  methods  yield  very 
similarly  shaped  functions.  The  main  difference  lies  in  the  scale  factor 
for  the  frequency  axis.  This  scale  factor  is  inversely  related  to  the 
effective  multipath  spread.  This  topic  will  subsequently  be  explored  in 
more  detail  to  obtain  a  reliable  estimate  of  the  multipath  spread  as  a 
function  of  path  parameters,  terrain,  and  weather. 
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C.  REFRACTIVE  INDEX  OF  THE  ATMOSPHERE 


The  density  of  the  atmosphere  and  its  radio  refractive  index  normally 
decrease  with  increasing  height.  A  generally  accepted  mathematical  model 
for  the  average  radio  refractive  index  as  a  function  of  height  is: 

N  =  Ns  e"ah 


where  Ns  =  surface  refractive  index 

a  =  scale  height  of  the  troposphere 
h  =  height  above  ground. 


The  agreement  of  this  model  with  various  sets  of  experimental  data  has  been 
investigated. 

One  set  of  experimental  data  consists  of  radiosonde  measurements  per¬ 
formed  in  San  Juan,  Puerto  Rico10  in  December  1965.  Measurements  made  at 
7:00  a. ir.  EST  on  7  December  1965  (Figure  5)  represent  typical  morning  con¬ 
ditions,  and  measurements  made  at  7:00  p.m.  EST  on  9  December  1965  (Figure 
6)  represent  typical  evening  conditions.  By  matching  the  slope  of  the 
experimental  curve  at  h  =  1  km,  it  can  be  seen  that 


N 


377  e 


-0 . 18h 


(16) 


gives  a  good  fit  to  the  experimental  data  over  the  primary  range  of  interest. 

The  corresponding  effective  earth  radius  can  be  found  from  the  equa¬ 
tion^-1  (also  see  Section  II-E) 


(17) 


where  Rq  is  the  actual  earth  radius  (6371  km).  For  N  given  by  equation  16, 
Re  =  9902  km. 


The  gradient  in  the  radio  refractive  index  near  the  earth's  surface 
can  also  be  found  from  the  formula  given  by  ESSA1^  based  on  the  empirical 
da  ta : 


AN 

Ah 


_  __  0.005577N 

1 .32  e  s 


(18) 
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For  Ng  =  377,  the  following  is  obtained: 


AN 

Ah 


59.9/km. 


When  the  corresponding  quantity  is  evaluated  directly  from  the  data  pre¬ 
sented  in  Figure  5,  the  following  is  obtained: 


AN 

Ah 


56/km. 


100  200  100  400 

IV  LVnts 


Figure  5.  Rawinsonde  Profile, 
San  Juan,  7  December  1965,  1200Z 
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Figure  6.  Rawinsonde  Profile, 
San  Juan,  9  December  1965,  0000Z 


The  agreement  between  the  experimental  data  and  the  results  based  on  equa¬ 
tion  18  is  very  good. 

For  the  typical  evening  situation  shown  in  Figure  6,  the  gradient  and 
the  value  of  the  radio  refractive  index  at  1  km  can  be  matched  to  obtain: 


N 


372  e 


-0.15 


h 


This  leads  to  an  effective  earth  radius  of  9176  km. 


(19) 


An  example  on  a  different  type  of  evening  data  is  shown  by  using  the 
data  obtained  from  7:00  p.m.  EST  on  13  December  1965  (Figure  7).  The  de¬ 
tailed  profile  is  quite  different  from  the  data  obtained  at  the  same  time 
of  day  on  9  December  1965.  Nevertheless,  a  good  fit  is  obtained  for: 


N 


-0.155  h 

362  e 


(20) 
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Figure  7.  Rawinsonde  Profile, 
San  Juan,  13  December 
1965,  0000Z 


The  corresponding  effective  earth  radius  is  approximately  the  same  as 
the  one  obtained  from  the  9  December  1965  data. 

According  to  typical  graphs,  the  effective  earth  radius  is  about  8 
percent  greater  in  the  morning  than  in  the  evening  for  this  over  water  path 
in  December.  Much  more  extreme  cases  can  be  found,  as  well  as  cases  where 
Re  is  smaller  in  the  morning  than  in  the  evening.  The  results  presented 
do,  however,  appear  to  be  typical. 

In  the  remaining  part  of  this  section  it  will  be  assumed  that  equation 
18  adequately  represents  the  gradient  of  the  refractive  index. 

D.  RAY  TRACING  IN  THE  ATMOSPHERE 

The  transmitted  signal  is  radiated  from  a  narrow  beam  antenna  in  a 
direction  more  or  less  tangential  to  the  earth's  surface.  The  density 
of  the  atmosphere,  and  therefore  its  radio  refractive  index,  decreases  with 
height.  This  decrease  in  radio  refractive  index  will  cause  the  transmitted 
beam  to  be  bent  toward  the  earth's  surface,  rather  than  to  travel  in  a 
straight  line.  This  bending  will  cause  a  change  in  the  sxze  and  location 
of  the  common  volume  of  the  intersection  between  the  transmitter  and  re¬ 
ceiver  beams.  It  is  the  purpose  of  this  section  to  evaluate  the  actual  ray 
paths  as  a  function  of  antenna  elevation  angle  and  surface  refractive  index 
for  a  commonly  used  model  of  tropospheric  density  variations. 

The  geometry  of  the  situation  is  shown  in  Figure  8.  The  beam  is 
launched  from  T  at  an  angle  y0  with  the  local  tangent  to  the  earth's  sur¬ 
face.  At  the  angle  /9  the  ray  is  h  meters  above  the  earth's  surface.  To 
solve  this  problem,  an  equation  for  ray  bending  will  be  set  up. 
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It  is  assumed  that  the  radio  refractive  index  of  the  atmosphere  de¬ 
creases  exponentially  with  height: 


n  =  1  +  n0  e  ar  .  (21) 

The  bending  of  the  ray  due  to  passage  through  an  infinitesimal  radial 
distance,  dr,  is  shown  in  Figure  9.  Using  Snell's  law  for  a  spherically 
symmetric  geometryl3f 

r  n  cos  y  =  constant  along  ray  =  k.  (22) 


Figure  9.  Differential 
Refraction  Geometry 


dr 

* 


Furthermore 


dr  ft  y  rd0  (23) 

may  be  derived  from  the  geometry.  By  use  of  equaticns  (21)  ,  (22)  ,  and 
(23) ,  it  follows  that 

dy  _  d-y  dn  dr 
dfl  -  dn  dr  d/3 

§£  =  1  -  ran  e“ar  .  (24) 

d/3  o 

This  equation  is  valid  for  y  <  10°,  a  condition  that  is  always  met  in  a 
practical  situation. 

Since  r  changes  only  by,  at  most,  a  few  km  in  several  thousand,  the 
right  hand  side  of  equation  (24)  is  assumed  to  be  constant.  We  thus  obtain 
the  equations : 


*2  = 

d/3 


K. 


(25) 


and 


_aD 

K_  =  1-Rane  o  (26) 

3  o  o 

_aR 

where  R  is  the  earth's  radius  and  -an  e  o  is  the  rate  of  change  of  the 
radio  refractive  index  near  the  earth's  surface.  When  equations  (23)  and 
(25)  are  combined,  equation  (27)  for  the  ray  path  is  obtained: 

r  =  Rq  exp  Tyo  P  +  0.5  K3p\  p  <5°  (27) 


The  distance  between  the  ray  and  the  earth's  surface  (for  a  smooth  earth) 
is  therefore : 


h  = 


R 

o 


+  0.5  K 

3 


P<  5° 


(28) 


To  use  equation  (28)  values  for  must  be  obtained.  From  equation  (26)  it 
can  be  seen  that: 


K  =  1  -  R  ^ 
3  o  Ar 


10 


-6 


(29) 


where  N  is  in  the  conventional  N  units. 


AN/Ar  can  be  found  by  actual  measurements  of  N  at  various  heights.  As  this 
is  not  practical  in  most  cases,  use  is  made  of  the  equation  given  by  ESSA^-2. 


AN  _  „  +0.005577N 

r —  =  -7.32  e  s 

Ar 


(30) 


where  N  is  the  surface  radio  refractive  index, 
s 


The  surface  refractive  index,  Ns ,  changes  from  hour  to  hour  during  the 
day  and  also  from  day  to  day  throughout  the  year.  As  the  surface  refractive 
index  changes  so  will  the  ray  path  and  the  location  and  size  of  the  common 
volume.  Ns  can  be  monitored  continuously,  but  in  most  cases  it  is  more 
practical  to  calculate  it  from  the  equation^ ; 

K  /  K  H  P  \ 

Ns  =37TTT\P+  273-TT/  (31) 

where 

Kx  =  77.6 
K2  =  48.1 

t  =  temperature  in  °C 
P  =  atmospheric  pressure  in  millibars 
-  H  =  relative  humidity  in  percent 

=  saturated  water  vapor  pressure  in  millibars 


The  only  uncommon  quantity  in  equation  (31)  is  P  .  Tables  of  P 
ily  available.  Over  a  limited  temperature  range  centered  on  15®C 
(32)  gives  quite  accurate  values  for  Pw: 


are  read- 
,  equation 


P 

w 


13. 2e 


0. 0391t 


millibars . 


(32) 


In  the  absence  of  anything  but  the  usual  meteorological  measurements 
for  the  location  and  time  in  question,  one  is  able  to  estimate  the  actual 
ray  path  for  a  proposed  troposcatter  circuit.  This  estimate  of  the  actual 
ray  path  will  enable  one  to  form  a  meaningful  estimate  of  the  size  of  the 
common  volume  and  thus  of  the  correlation  bandwidth. 

It  is  assumed  that  beams  are  launched  at  various  elevation  angles  at  the 
transmitting  site.  The  ray  path  is  then  computed  based  on  the  above  theory. 

The  vertical  size  of  the  common  volume  due  to  an  antenna  of  a  given  effective 
beamwidth  is  the  difference  in  height  of  the  two  ray  paths  at  the  midpath  point. 
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Tables  that  can  be  used  to  predict  the  size  of  the  common  volume  and 

corresponding  multipath  spread  for  large  combinations  of  take-off  angles , 

path  length,  and  surface  refractive  index  have  been  presented  in  the  first 

interim  report  (Ref.  6).  Based  on  these  tables,  plots  of  ray  path  versus 

beta  in  degrees  are  presented  in  Figures  10  through  13  to  illustrate  the 

relative  refraction  of  the  rays  for  different  values  of  surface  refractive 

index,  N  . 

s 

This  ray  tracing  technique  has  been  used  to  computer  calculate  the 
multipath  spread  given  the  local  weather  conditions  (to  find  Ns) ,  the  path 
distance,  the  two  take-off  angles  based  on  4/3  earth  radius,  and  the 
effective  beamwidth  of  the  antennas.  The  concept  of  effective  beamwidth 
is  discussed  in  subsection  H.  The  computer  program  is  listed  in  Appendix  C 


Initial  Take  Off  Angle  Gamma  *  0  Degree 


Ns  •  301 
Ns  ■  350 


Beta  !.n  Degrees 


Figure  10.  Height  of  Ray  versus  Beta  in  Degree. 


Height  of  Ray  Above  Ground  -  Meters 


Height  of  Ray  Above  Ground  -  Meters 


Arc  Length  in  Kilometers 

Figure  13.  Height  of  Ray  versus  Arc  Length  in  km.  Initial  Take-Off  Angle 

Gamma  =  1.0  Degree 
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E.  EFFECTIVE  EARTH  RADIUS 

As  shown  in  the  previous  section,  the  transmitter  rays  are  bent  toward 
the  earth  due  to  the  decrease  in  radio  refractive  index  with  height.  Ordi¬ 
nary  geometric  ray  tracing  can,  therefore,  not  be  directly  used  to  predict 
the  location  and  size  of  the  common  volume. 

An  approximate  method  for  circumventing  this  difficulty  is  to  assume 
that  the  rays  travel  in  straight  lines  but  that  the  earth  has  an  effective 
radius,  Re,  greater  than  the  actual  one.  This  effective  radius  is  chosen 
so  that  the  height  of  the  common  volume  above  ground  level,  when  the  rays 
are  assumed  straight,  is  the  same  as  it  would  be  in  the  actual  case.  The 
geometry  of  the  situation  is  explained  in  Figure  14.  The  geometrical 
angular  distance  between  the  midpath  point  and  the  transmitter  is  /8 .  The 
effective  angular  distance  is  £  .  When  yQ  *  0,  from  equation  (28),  h  becomes 

o  cif 

h  (  y  *  0)  =  R  (e  u * 3  3  M  -1)  .  (33) 

o  o 


This  equation  for  the  fictitious  earth  will  now  be  modified, 
straight  line  ray  corresponds  to  K  ■  1,  it  can  be  seen  that 


h  = 


R 

e 


-1). 


Since  a 


(34) 


Real  Earth  with  Tropospheric 
Refraction 


Fictitious  Earth  with 
Larger  Radius 
but  No  Refraction 


Figure  14.  Geometry  Showing  Relation  Between  Real  and  Effective  Earth  Radius 
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Since  the  linear  distance  between  the  transmitter  and  the  midpoint  of  the 
path  is  unaffected  by  atmospheric  conditions. 


(35) 


Thus 


R 

o 


,  0.5K, 
(e  3 


02_ 


1) 


R 

e 


(36) 


A  numerical  solution  to  equation  (36)  can  be  obtained  by  use  of  a  computer. 
However,  since  P  <  5°  and  K3  and  R0/Re  are  less  rhan  one  we  can  approximate 
the  exponential  function  by  the  first  two  terms  of  its  power  series  expan¬ 
sion.  When  this  is  accomplished  equation  (37)  is  obtained: 


(37) 


When  K3  is  expressed  in  terms  of  its  constituents,  the  final  expression  for 
the  effective  earth  radius  in  terms  of  the  parameters  of  the  atmosphere  is 
obtained : 


R 

e 


(38) 


F.  RAKE  MEASUREMENT 

One  method  for  studying  the  structure  of  the  common  volume  in  detail 
is  the  RAKE  technique.  This  method  is  similar  to  the  method  used  in  radar. 

A  signal  that  can  be  resolved  well  in  both  time  and  frequency  is  transmit¬ 
ted.  The  received  signal  is  resolved.  Differences  in  time  shift  of  various 
components  of  the  received  signal  correspond  to  differences  in  propagation 
delay.  Differences  in  frequency  correspond  to  doppler  shift  suffered  by 
the  signal  during  transmission. 

Results  obtained  in  February  1965  on  the  484  km  path  from  Fort  Bragg, ^ 
North  Carolina  to  a  site  located  about  20  km  Northeast  of  Washington,  D.C. 
will  be  studied  now.  The  most  important  parameters  of  this  path  are  listed 
below.  It  is  assumed  that  the  signals  travel  in  straight  lines  and  that  the 
effective  earth  radius  is  4/3  of  the  actual  one  obtained  as  shown  in  Figure 
15.  This  Figure  shows  the  earth  surface,  the  grazing  ray,  and  the  family 
of  confocal  ellipsoids  corresponding  to  various  time  delays  relative  to  the 
transmission  delay  of  the  grazing  ray. 
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Operating  frequency 
Transmitter  power 
Path  length 


Summary  of  Some  Path  Parameters 


Antenna  Size 


910  MHz  Shell  thickness  corres-  200  m 

ponding  to  single  tap 
8  kW  near  top  of  common  volume 

484  km  Shell  thickness  corres-  550  m 

ponding  to  single  tap  near 
8.5  m  bottom  of  common  volume 


Beamwidth  3  deg 

Height  of  bottom  of  3.4  km 

scattering  volume 

Height  of  top  of  16  km 

scattering  volume 


Differential  time  delay  4.6  ps 

between  top  and  bottom  of 
common  volume 

Effective  scattering  3.26  deg 

angle  (4/3  earth) 


Figure  15.  Longitudinal  and  Transverse  Sections  of  Ellipsoids  of 

Constant  Time  Delay 


The  antennas  were  aligned  with  beam  intersection  as  shown  by  the  cir¬ 
cle  in  Figure  15.  The  scattering  function  displays  the  spreading  of  the 
transmitted  signal  in  time  and  frequency  as  shown  in  Figure  16.  The  fre¬ 
quency  spreading  is  due  to  doppler  shift  due  to  crosspath  wind.  This  dop- 
pler  shift  affects  the  fading  rate  but  not  the  frequency  correlation  func¬ 
tion.  Therefore,  it  will  not  be  pursued  further  except  where  it  can  be 
used  as  an  aid  in  locating  the  effective  common  volume  by  a  comparison  of 
doppler  shifts  and  crosspath  wind  data.  The  time  spreading  is  of  primary 
interest.  Since  the  vertical  axis  through  the  antenna  pattern  in  Figure  15 
is  labeled  in  relative  time  delay,  a  direct  mapping  from  the  observed  scat¬ 
tering  function  to  the  common  volume  can  be  performed. 

Due  to  the  small  number  of  records  available,  it  is  difficult  to  draw 
definite  conclusions  about  the  nature  of  the  observed  scattering  functions. 
However,  if  this  limitation  is  kept  in  mind  some  tentative  conclusions  may 
be  drawn. 

The  scattering  function  shown  in  Figure  16  was  measured  at  11:57  a.m. 
on  17  February  1965;  the  one  in  Figure  17  at  1:37  p.m.  on  19  February  1965; 
and  the  one  in  Figure  18  at  5:26  p.m.  on  16  February  1965.  If  one  assumes 
that  the  prevailing  weather  conditions  were  similar  on  the  days  of  the 
measurements,  the  differences  in  scattering  functions  must  be  related  to 
the  differences  in  time  of  day  when  the  measurements  were  performed.  If 
these  differences  are  typical,  the  multipath  spread  tends  to  be  relatively 
narrow  in  the  morning  and  increases  during  the  day. 


( 


Frequency  -  Hz 

Figure  16.  Scattering  Function,  Record  No.  142,  11:57  a.m., 17  February  1965 
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If  it  is  assumed  that  the  bending  of  the  transmitted  and  received  rays 
can  be  accounted  for  by  a  4/3  earth  radius  throughout  the  day,  the  multi- 
path  spread  on  the  antenna  pattern  can  be  plotted  as  shown  in  Figures  19, 

20,  and  21  for  various  times  of  the  day.  Significantly,  it  can  be  seen  that 
at  least  for  this  fairly  long  path,  the  effective  common  volume  as  measured 
by  this  RAKE  technique  is  always  much  smaller  then  the  one  predicted  by  the 
intersection  of  the  transmitter  and  receiver  antenna  beams. 

A  point  still  open  to  question  is  the  height  of  the  effective  common 
volume  above  ground  level.  This  point  can  be  resolved  if  there  is  a  rela- 
strong  known  crosspath  wind  at  the  common  volume.  Birkemeier^ 
reports  on  measurements  made  on  a  230  km  path  between  Cedar  Rapids,  Iowa 
and  Arlington,  Wisconsin.  Birkemeier  measured  the  doppler  shift  on  the 
Continuous  Wave  signal  as  a  function  of  antenna  azimuth.  At  the  same  time 
he  observed  the  crosspath  wind.  During  one  measurement,  the  wind  was  South¬ 
easterly  below  and  Northwesterly  above  1600  meters.  The  3  dB  antenna  beam- 
width  corresponded  to  a  common  volume  height  range  from  800  to  3700  meters 
above  ground  level,  based  on  a  4/3  earth  radius.  Birkemeier  found  that 
the  sign  of  the  average  doppler  shift  corresponded  to  the  direction  of  the 
wind  below  1600  meters.  Thus,  the  lower  800  meters  of  the  common  volume 
transferred  the  majority  of  the  transmitted  signal  power. 


Figure 

Time 


19.  Longitudinal  and  Transverse  Sections  of  Ellipsoids  of  Constant 
Delay  Showing  Effective  Common  Volume,  11:57  a.m.,  17  February  1965 
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Figure  20.  Longitudinal  and  Transverse  Sections  of  Ellipsoids  of  Constant 
Time  Delay  Showing  Effective  Common  Volume,  1:37  p.m.,  19  February  1965 


Time  Delay  Showing  Effective  Common  Volume,  5:26  p.m.,  16  February  1965 
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In  an  experiment  on  25  July  1966,  the  average  doppler  shift  on  the 
received  signal  corresponded  to  a  crosspath  wind  of  10  meters  per  second. 
Figure  22  shows  the  measured  wind  profile  versus  height  for  the  morning  of 
25  July  1966.  The  wind  velocity  was  in  the  10  meter  per  second  range  from 
about  1200  to  1500  meters  above  ground  level  indicating  that  this  is  the 
location  of  the  effective  common  volume. 

These  two  detailed  investigations  tend  to  reinforce  the  contention 
that  the  scattering  angle,  rather  than  antenna  beamwidth,  usually  tends  to 
restrict  the  size  of  the  common  volume.  Furthermore,  from  Figure  21,  it  is 
apparent  that  the  effective  antenna  beamwidth  at  least  for  this  occasion 
was  only  one-half  of  the  actual  beamwidth  with  the  majority  of  forward 
scatter  occurring  in  the  lower  portion  of  the  common  volume  formed  by  the 
intersection  of  the  transmit-receive  beams. 


20  16  12  8  4  0  4  8 

Crosspath  Wind  Speed  (m/sec) 

Figure  22.  Tropospheric  Wind  Speed  Profile 


G.  STABILITY  OF  THE  AIR  IN  THE  COMMON  VOLUME 


The  stability  of  the  air  within  the  common  volume  plays  an  important 
part  in  determining  the  correlation  bandwidth  of  the  troposcatter  link. 
Frequently  under  unstable  conditions  the  correlation  bandwidth  has  been 
found  to  drop  by  a  factor  of  approximately  two.  It  is,  therefore,  import¬ 
ant  to  ascertain  the  meteorological  parameters  contributing  to  instability. 

The  reason  for  instability  of  the  air  in  the  common  volume  must  be 
sought  in  the  thermodynamics  of  the  air  in  the  layer  from  about  500  to 
2000  meters  above  smooth  earth  level.  When  a  small  volume  of  air  is 
raised  it  will  undergo  an  adiabatic  expansion.  Work  is  required  to  make 
the  air  expand  adiabatically . 
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As  the  air  expands  it  also  cools.  If  it  were  an  ideal  gas,  this 
could  go  on  indefinitely.  However,  since  the  air  contains  water  vapor  it 
will  eventually  reach  a  point  where  the  temperature  equals  the  dew  point 
of  the  mixture.  Further  elevation  of  the  air  will  cause  an  adiabatic 
expansion  with  condensation,  a  so  called  pseudoadiabatic  expansion.  In 
this  latter  phase,  energy  is  released. 

If,  on  balance,  energy  is  released  when  a  small  volume  of  air  is 
raised,  the  air  is  in  an  unstable  condition. 

Estimation  of  the  potential  stability  of  the  air  can  conveniently  be 
done  by  use  of  curves  of  temperature  and  dew  point  as  a  function  of  pres¬ 
sure.  Such  charts  are  prepared  twice  daily  by  a  number  of  weather  fore¬ 
casting  stations  throughout  the  world.  We  will  give  an  example  of  how 
such  a  chart  can  be  used  to  estimate  the  potential  stability  of  the  air 
including  the  common  volume. 

Let  us  consider  the  pressure  temperature  curve  shows  in  Figure  23. 

We  will  estimate  the  work  required  to  lift  a  small  volume  of  air  from  the 
950  mb  level  where  it  will  be  at  a  temperature  to  the  850  mb  level 
where  it  will  be  at  the  temperature  T^. 

Initially  the  air  will  expand  adiabatically  parallel  to  the  adiabatic 
direction  until  it  is  saturated.  This  is  the  point  M2.  From  M2  to  the 
850  mb  level  the  air  will  expand  pseudoadiabatically ,  i.e.,  with  condensa¬ 
tion.  The  air  thus  will  arrive  at  the  850  mb  level  at  the  point  M3 .  The 


28 


area  between  the  P— T  curve  and  the  path  M^,  M 2 ,  M3,  and  M^,  counting  area 
in  the  right  side  of  the  curve  as  positive,  represents  the  work  required 
to  move  a  small  volume  of  air  from  the  950  mb  level  to  the  850  mb  level. 

If  this  work  is  positive,  the  atmosphere  is  stable,  if  it  is  negative 
the  atmosphere  is  unstable. 

This  method  allows  estimation  of  the  stability  of  the  air  in  the 
common  volume  when  sonde  data  is  available.  If  one  postulates  a  suitable 
thermodynamic  model  for  the  atmosphere,  it  is  possible  to  obtain  a  reason¬ 
able  estimate  of  the  P— T  curve  from  a  sequence  of  ground  based  temperature 
measurements  even  in  the  absence  of  sonde  data . 


In  the  following  we  will  assume  that  there  is  no  inversion  layer 
between  the  ground  and  the  comnon  volume.  In  the  summer  the  air  above  the 
ground  will  heat  during  the  day  due  to  conduction  and  radiation  from  the 
earth.  This  warm  air  will  gradually  rise  and  supply  the  atmosphere  with 
warmer,  and  frequently  more  moist  air.  If  we  assume  that  the  P-T  curve 
was  vertical  in  the  morning  it  will  acquire  a  decreasing  negative  slope 
throughout  the  day.  The  slope  will  be  least  at  the  time  of  day  when  the 
ground  level  temperature  is  highest.  This  situation  is  shown  in  Figure 
24. 


We  see  from  Figure  24  that  if  the  P-T  curve  is  as  assumed,  an  unstable 
condition  exists  in  the  common  volume. 


Assumed 
P-T  Curve 


Figure  24.  P-T  Curve  Showing  Atmospheric  Instability 
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From  this  we  can  conclude  that  under  normal  circumstances  the  greatest 
chance  for  turbulence  occurs  at  about  the  time  when  the  ground  level  temp¬ 
erature  reaches  it  maximum  value. 

Since  the  degree  of  turbulence  has  a  direct  effect  on  multipath  spread, 
a  scale  factor  can  be  associated  with  the  stability  of  the  troposphere  in 
and  near  the  coitmon  volume.  This  scale  factor  will  vary  from  unity  in  the 
case  of  complete  turbulence  to  as  small  as  0.4  in  the  case  of  normal 
stratification.  Under  unusual  propagation  conditions,  such  as  are  some¬ 
times  found  over  water,  the  scale  factor  can  attain  even  smaller  values. 

The  probability  distribution  of  the  stability  of  the  troposphere  is 
therefore  directly  related  to  the  probability  distribution  of  the  multi- 
path  speed  and  in  turn  directly  relatable  to  the  correlation  bandwidth  of 
the  troposcatter  channel. 

H.  DEVELOPMENT  OF  THE  COREELATION  BANDWIDTH  MODEL 


Within  the  framework  of  the  foregoing  discussion,  an  empirical  model 
has  been  developed  which  successfully  predicts  the  envelope  frequency 
correlation  function  for  the  series  of  field  tests  that  were  made  during 
this  program.  It  is  believed  that  this  model  has  wider  applicability  in 
that: 

1_  There  is  a  plausible  theoretical  hypothesis  to  support  it,  and 

2_  Application  to  paths  not  measured  here  has  yielded  good  agreement. 

The  approach  taken  has  been  to  obtain  a  model  for  the  lower  decile  of 
the  measured  correlation  functions,  that  is,  those  with  maximum  multipath 
spread.  A  study  of  the  correlation  function  plots  in  a  following  section 
has  revealed  that  to  a  good  approximation,  the  upper  decile  of  measured 
correlation  functions  corresponds  to  those  of  the  lower  decile  with  a 
frequency  scale  factor  of  2.5.  Therefore,  emphasis  has  been  placed  on  the 
prediction  of  the  lower  deciles.  It  is  not  crucial  to  this  model  but  it 
has  been  argued  that  the  lower  deciles  correspond  to  those  occasions  when 
turbulent  scattering  has  been  predominant,  the  upper  deciles  occurring 
during  conditions  of  stratification.  The  preceding  subsection  has  detailed 
how  a  prediction  of  stability  or  turbulence  in  the  common  volume  is  possible 
with  the  general  conclusion  that  more  turbulence  is  to  be  expected  when  the 
surface  temperatures  are  maximal  as  would  normally  occur  during  the  after¬ 
noon.  Verification  of  this  hypothesis  has  been  made  for  a  limited  number 
of  individually  measured  tests. 

The  correlation  function  has  been  modeled  as  Gaussian.  As  just  dis¬ 
cussed,  any  of  the  several  earlier  models  indicate  that  this  is  a  valid 
approximation  if  a  suitable  scale  factor  (or  equivalently  variance)  is 
found.  Moreover,  the  results  of  the  tests  included  here  point  up  that  in 
actuality,  the  correlation  function  is  nearly  Gaussian  in  the  region  of 

interest,  0.4  <  p  <1. 

*e 
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It  was  found  quite  early  that  the  actual  half  power  antenna  beamwidth 
(HPBW)  had  little  bearing  on  the  measured  correlation  functions.  If  it 
had,  the  X-  and  C-band  tests  with  differing  antenna  beamwidths  would  have 
exhibited  considerably  different  correlation  bandwidths,  while  in  fact 
the  deciles  have  been  found  nearly  equivalent.  This  same  conclusion  with 
measurements  at  a  single  frequency  but  differing  antenna  sizes  has  been 
demonstrated  previously 7. 

It  appears  that  in  most  practical  links  of  path  length  over  100  miles, 
the  intersection  of  the  transmit  and  receive  beam§  is  considerably  differ¬ 
ent  than  that  part  of  the  atmosphere  which  efficiently  produces  forward 
scatter.  Consequently,  in  the  calculation  of  the  multipath  spread,  it  is 
convenient  to  introduce  the  concept  of  effective  beamwidth  «e  which  is 
the  angular  extent  of  the  beam  producing  effective  scattering.  Clearly 
the  effective  beamwidth  should  be  limited  to  approximately  the  HPBW  of 
either  antenna  a .  It  is  assumed  that  a  is  the  same  for  both  the  transmit 
and  receive  antennas  as  it  was  for  these  tests.  Since  a  will  be  the  key 
parameter,  this  is  not  a  critical  assumption.  In  an  earlier  reference  to 
a  paper  by  Gordon^,  it  was  mentioned  that  the  effective  beamwidth  is  often 
taken  to  be  the  smaller  of  the  actual  beamwidth  or  one-half  the  scatter 
angle.  In  fact,  in  the  three  earlier  interim  reports  (References  6,  15, 
and  20)  this  assumption  was  made.  In  order  to  fit  the  data,  this  un¬ 
fortunately  led  to  a  "variable"  scale  factor  in  the  exponential  model  for 
each  of  the  four  paths.  The  model  developed  here  is  essentially  the  same 
as  given  in  the  earlier  interim  reports  but  it  is  assumed  more  generally 
that  the  effective  beamwidth  is  equal  to  the  product  K(6)  •  0  where  K(6) 

has  been  empirically  determined.  Thus  the  variation  between  the  four  paths 
will  be  taken  out  of  the  scale  factor  (Sfc)  and  be  placed  more  naturally  in 
the  effective  beamwidth. 

Toward  the  development  of  this  model,  the  envelope  correlation  function 
can  be  represented  as 


Pe  =  exp  [  -  (st  A  f)2]  (39) 

where  S^-  is  a  scale  factor  for  turbulent  scattering  to  be  determined,  A  is 
the  multipath  spread  as  calculated  by  the  ray  trace  model  for  an  effective 
beamwidth  ae,  and  f  is  the  frequency  separation.  Ss,  the  scale  factor  for 
a  stable  common  volume,  will  be  S  /2.5  or  0.4  S  . 

For  a  given  multipath  spread,  geometrical  considerations  imply  that 
the  effective  beamwidth  must  increase  for  shorter  paths.  Since  the 
observed  multipath  spreads  were  not  markedly  different  over  the  four  paths, 
it  follows  that  the  effective  beamwidth  must  be  greatest  for  the  shortest 
of  the  paths,  that  is.  Model  City  -  Ontario  Center.  To  insure  that  the 
effective  beamwidth  ae  was  never  greater  than  the  actual  value  of  a  ,  the 
value  of  ae  for  this  path  was  set  equal  to  the  smallest  actual  value  of  or, 
the  X-band  figure  of  1.07  degrees.  It  is  realized  that  this  assumption  is 
somewhat  arbitrary  since  even  for  this  short  path  (85  miles)  there  is  no 
firm  assurance  that  the  entire  beamwidth  produces  scatter.  It  is,  however, 
a  reasonable  working  hypothesis  and  should  be  non-crucial  for  any  longer 
paths  (or  more  precisely,  larger  scatter  angles) . 
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With  the  assumed  value  of  ae  for  the  Ontario  Center  link  it  is 
possible  to  find  the  multipath  spread.  A  value  of  Ns  must  be  adopted  and 
in  the  following,  unless  otherwise  stated,  Ns=  301  has  been  used  exclus¬ 
ively  implying  the  usual  4/3  earth  radius  approximation.  The  validity  of 
this  approach  will  be  discussed  later. 

The  calculation  of  the  multipath  spread  A  was  here  accomplished  using 
the  ray  trace  program  listed  in  Appendix  C.  Other  than  Ng  and  a  .  A  depends 
on  the  path  distance  d,  the  two  take-off  angles,  y £  and  yw,  an§  the 
altitude  of  the  two  terminals.  These  are  available  from  the  contour 
profiles  shown  in  Section  III.  The  notation  is  shown  in  Figure  25.  The 
calculated  value  of  A  for  the  Ontario  Center  link  is  found  to  be  0.157 
microseconds.  To  agree  with  the  X-band,  Ontario  Center,  summer  lower 
decile  correlation  bandwidth  of  1.5  MHz  (at  P&  =  0.4),  St  must  then  be 
equal  to  4.06. 


St  =  4.06 

and  therefore  Ss  *  1-62  • 


Figure  25.  Geometry  for  Empirical  Troposcatter  Model 


It  is  of  interest  that  this  value  of  St  is  not  radically  different  from 
Rice's  factor  of  2  tt.  Additionally,  one  series  of  simulator  measurements, 
in  which  A  was  progranmed  in  with  the  correlation  bandwidth  being  measured , 
indicated  an  St  of  about  2.5  (Reference  21,  p  258). 

If  St  is  assumed  fixed  for  the  three  remaining  paths,  it  is 
possible  to  find  the  multipath  spread  necessary  to  match  the  measured 
lower  decile  correlation  bandwidths.  These  multipath  spreads,  in  turn, 
through  the  ray  trace  program,  imply  an  effective  beamwidth  necessary  to 
produce  than. 

Thus  for  the  desired  value  of  St  a  determination  of  the  value  of 
effective  beamwidth  can  be  made  for  each.  These  values  are  plotted  versus 
the  scatter  angle, 


(40) 

(41) 


6  =  YE  +  YW 


(42) 


in  Figure  26.  These  scatter  angles  are  based  on  Ng  =  301  (Re  =  £  Rq)  . 
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i. 


0  Scatter  Angle  In  Degrees 


Figure  26.  Empirical  Effective 
Beamwidth  for  Application  in  Ray 
Tracing  Computer  Program 


The  analytic  function 

a  =  ( 0.57(9  -  0.77)"0-45  -  0.36]  02  ( a  ,  0  in  degrees)  (43) 

e  i  j  e 

gives  a  reasonable  fit  to  the  experimental  points.  Equivalent  to  this 
figure.  Figure  27  gives  K(0)  as  a  function  of  the  scatter  angle  0,  where 
K (0) •  0  is  the  effective  beamwidth.  That  is, 

K (0)  =  0.57  0  (0  -  0.77)“0,45  -  0.36  0  (44) 


Figure  27.  k  as  a  Function  of 
Scatter  Angle,  where  K  is  the 
Quotient  of  the  Effective  Beam- 
width  and  the  Scatter  Angle 


o.; 


These  figures  illustrate  that  in  general,  if  the  scale  factor  St  is 
assumed  constant,  cte  is  not  equal  to  one-half  the  scatter  angle;  that  is, 

K  #  0.5  in  general .  It  is  also  apparent  that  a  e  may  become  larger  than 
a,  the  actual  half -power  beamwidth.  In  practice,  &e  must  be  limited  to 
a. 

The  scatter  angle  9  is  felt  to  be  the  single  most  important  parameter 
in  characterizing  a  link  since  it  depends  on  both  the  path  distance  and 
the  take-off  angles  {as  given  by  equation  42).  It  is  likely  that  the 
effective  beamwidth,  and  ultimately  the  correlation  bandwidth,  depends 
independently  on  the  path  distance  and  take-off  angles  in  a  way  that  is 
more  complicated  than  summarized  in  the  single  parameter,  scatter  angle. 
This  could  be  an  area  of  further  research.  Nevertheless,  lacking  further 
data,  there  is  little  question  that  a  fairly  accurate  model  of  a  path  can 
be  made  from  the  single  concept  of  an  effective  beamwidth  which  is  a 
function  only  of  the  scatter  angle. 

I 

The  calculation  of  the  scatter  angle  has  been  predicated  on  the 
assumption  of  Ns  =  :>01  since  this  provides  a  base  for  future  predictions. 
This  implies  that  the  scatter  angle  is  to  be  considered  a  fixed  parameter 
for  a  specific  link,  Tr.  actuality  this  condition  will  rarely  be  realized; 
for  the  links  tested  here  Ng  varied  from  approximately  301  to  370.  Never¬ 
theless,  for  predicting  the  correlation  bandwidth  for  other  links,  the 
scatter  angle  based  on  4/3  earth  radius  (Ns  =  301)  should  also  be  used. 

Thus,  even  roo  .gh  it  is  very  convenient  to  characterize  a  link  by  a 
fixed  scatter  ar^gle  it  must  be  kept  in  mind  that  the  actual  value  of  0  is 
dependent  on  o.  Lime  varying  Re  (and  to  a  lesser  degree  of  importance  on 
time  varying  Ye,  Yw)  -  This  fact  will  lead  to  negligible  modifications  for 
the  case  of  longer  paths  where  9  is  large  but  it  is  to  be  noted  that  for 
short  paths  (small  9)  ,  K  is  very  sensitive  to  changes  in  9.  Something 
approaching  this  behavior  was  noted  for  the  Ontario  Center  path  where  the 
short  term  variations  (periods  of  1  to  3  hours)  were  more  pronounced  than 
for  the  other  three  longer  links.  Ideally,  figures  such  as  Figure  27 
could  be  made  for  a  set  of  Ns  values.  Here,  the  method  of  computer  re¬ 
duction  did  not  lend  itself  to  this  type  of  analysis  and  it  must  be  con¬ 
sidered  an  area  of  future  study. 

It  has  become  evident,  however,  that  even  in  the  region  of  Figure  27 
where  K  is  sensitive  to  9,  the  observed  variations  between  the  upper  and 
lower  decile  correlation  bandwidths  is  much  too  great  to  be  accounted  for 
by  simple  changes  in  Ng  (and  therefore  0  and  ae)  .  This  is  felt  to  be  an 
important  point  and  lends  further  substance  to  the  argument  that  between 
the  upper  and  lower  deciles  there  is  a  change  in  the  common  volume  more 
severe  than  could  be  explained  by  variations  in  Ns;  that  is,  it  supports 
the  theory  of  turbulent  versus  layer  scattering. 

Although  Ng  is  considered  constant  in  the  determination  of  0,  N  also 
enters  as  a  parameter  in  the  ray  trace  procedure .  Here  it  may  be  alfowed 
to  vary  and  set  to  any  expected  value.  However,  as  will  be  discussed  in 
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the  next  subsection,  changes  in  Ns  in  the  ray  trace  program  have  only  a 
small  effect  on  the  correlation  bandwidth. 

It  will  be  instructive  to  illustrate  the  procedure  to  be  used  in 
predicting  the  correlation  bandwidth  for  a  hypothetical  link.  It  is 
assumed  that  d,  Y E»  Yw»  and  Ng  (or  P,  T,  %H)  are  known.  The  true  beam- 
width  will  also  be  known  but  it  plays  no  role  unless  ae  is  found  to  be 
greater  than  a  in  which  case  a  would  be  used. 

First,  the  scatter  angle  9  would  be  calculated  by  equation  42  in 
which  Ye»  Yw»  and  d/Re  are  all  based  on  the  approximation  of  4/3  earth 
radius.  Using  Figure  26  or  equation  43  the  effective  beamwidth  would  be 
found  and  adopted  unless  it  were  larger  than  the  actual  HPBW.  The  value 
of  Ng  would  be  calculated  from  the  expected  values  of  pressure,  tempera¬ 
ture,  and  relative  humidity  by  equation  31.  Table  I  in  the  next  sub¬ 
section  gives  typical  values  of  Ns  for  various  combinations  of  these 
parameters.  At  this  point,  computer  techniques  for  ray  tracing  greatly 
facilitate  the  calculation  of  the  multipath  spread  though  the  calculation 
could  be  approximated  by  hand.  A  listing  of  the  ray  trace  program  based 
on  the  derivations  in  Section  II-D  is  given  in  Appendix  C.  Knowing  A  , 
the  model  correlation  function  is  given  by  equati  '9,  that  ir» , 

pe  =  exp[  -(4.06  A  f)2]  (45) 

This  pe  is  explicitly  for  the  lower  decile  of  expected  correlation  band- 
widths.  The  upper  decile  is  of  the  same  form,  but  with  St  replaced  by 
Ss  =  1.62. 


I.  EFFECT  OF  WEATHER  ON  CORRELATION  BANDWIDTH  MODEL 

The  surface  index  of  refraction  Ns  may  be  calculated  from  the  tem¬ 
perature,  relative  humidity,  and  barometric  pressure  measured  at  the 
receiver  or  transmitter  site  by  equation  31.  Table  I  displays  the  values 
of  Ns  calculated  for  various  typical  values  of  these  parameters.  Over  the 
range  of  weather  conditions  shown  in  this  table  Ns  ranges  from  304  to  354. 
During  the  test  period  Ns  values  as  high  as  370  were  calculated  from  actual 
weather  conditions. 

As  pointed  out  previously,  variations  in  Ng  produce  two  changes  in  the 
empirical  model.  The  first  is  the  dependence  on  N  in  the  ray  trace  pro¬ 
gram.  There  the  rays  will  follow  slightly  different  paths  for  varying  Ng. 
Secondly,  the  scatter  angle  will  also  in  actuality  depend  on  Ns  through  a 
dependence  on  Rg.  This  changing  scatter  angle  will  affect  the  effective 
beamwidth  and  thus  the  correlation  bandwidth.  In  the  model  presented  here 
predictions  were  based  on  a  4/3  earth  radius  scatter  angle  and  this  second 
effect  will  at  this  point  be  considered  only  in  a  qualitative  manner. 

To  determine  the  first  of  these  effects, the  correlation  bandwidth 
model  was  evaluated  setting  Ng  equal  to  301  and  370  in  the  ray  trace  pro¬ 
gram  for  both  the  shortest  path.  Model  City-Ontario  Center,  and  the 
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longest  path.  Model  City-Whitford  Field.  Figures  28  and  29  are  plots  of 
the  resulting  predicted  correlation  coefficients.  It  is  evident  that 
changes  of  Ng  in  the  ray  trace  program  have  little  effect  on  the  corre¬ 
lation  bandwidth . 


Table  I 

CHANGES  IN  INDEX  OF  REFRACTION,  N  ,  WITH  CHANGES 
IN  TEMPERATURE,  RELATIVE  HUMIDITY,  AND® BAROMETRIC  PRESSURE 


Temperature 
Degrees  F. 

Relative 

Humidity 

Percent 

Barometric 

Pressure 

Inches  Mercury 

N 

s 

15 

60 

29.75 

326.0 

45 

60 

29.75 

328.9 

75 

60 

29.75 

348.6 

45 

30 

29.75 

303.9 

45 

60 

29.  75 

328.9 

45 

90 

29.75 

353.9 

45 

60 

29.50 

326.5 

45 

60 

29.75 

328.9 

45 

60 

30.00 

333.2 

Figure  28.  Change  in  Correlation 
Function  with  Change  in  Ng  in  Ray 
Trace  Program.  Ontario  Center 
Receive  Site. 


Figure  29.  Change  in  Correlation 
Function  with  Change  in  Ns  in  Ray 
Trace  Program.  Whit ford  Field 
Receive  Site. 
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Secondly,  as  mentioned  in  Section  II-H,  the  scatter  angle  0  will  in 
fact  be  time  variant.  For  6  written  as 

0  =  YE  +  YW  + 

it  is  possible  to  show  the  effect  of  a  changing  Ns  on  each  of  the  three 
components.  Table  II  gives  the  values  of  yE  ,  yw,  and  d/Re  for  the  four 
test  paths  under  conditions  of  Ng  equal  to  301,  350,  and  370.  The  take-off 
angles  vary  with  Re  in  the  following  way : 


o  e 


where  hg  and  h  are  the  elevations  of  the  site  and  obstruction  respectively, 
and  dQ  is  the  distance  from  the  site  to  the  obstruction.  The  values  of  0 
and  ae  (both  in  degrees)  for  the  three  conditions  are  also  listed.  The 
variation  of  0  is  seen  to  be  small  and  will  lead  to  only  minor  modifi¬ 
cations  in  the  effective  beamwidth  and  the  predicted  correlation  bandwidth 
for  three  of  the  paths.  For  the  Ontario  Center  link  ae  will  remain  1.07 
degrees  since  smaller  values  of  0  than  0.843  imply  ae>  a,  a  condition  for 
which  the  model  requires  a  =  a.  For  the  Point  Petre  link  the  predicted 
change  in  ae  will  be  large  but  still  insufficient  to  explain  the  observed 
variation  between  the  upper  and  lower  deciles.  This  variation  is  shown 
in  Figure  30.  For  the  remaining  two  paths  the  variation  is  very  much 
smaller. 


Table  II 

VARIATION  IN  SCATTER  ANGLE  WITH  N 

s 

(all  angles  in  degrees) 


Ns  =  301 

Re  =  8493 

km 

Ns 

R 

=  350 

=  9487 

km 

N 

b 

=  370 

=  10, 

06M  km 

VE 

Yw  d/R 

e 

8 

a 

e 

V 

d/R 

n 

i- 

nt ji Center 

-J. 124 

0.025  0. 94  3 

0.843 

1.07 

-0. 121 

■j.u3o 

0. 844 

0.753 

i .  r. 

-0.120 

0.0  32 

u .  79( 

0.  7C*b 

1..7 

Point  Petre 

-0.059 

-0.034  1.103 

1.010 

0.7  35 

-o.  05t> 

-0. J3u 

0.  987 

0.  901 

0.854 

-0.055 

-e.‘»28 

0.  931 

0.  848 

1 . 06 

Port  Byron 

-0.01  ’ 

0.049  1.296 

1. 3Jr 

U.  680 

-U. 005 

u.o52 

1 .  lbu 

1.207 

0.682 

-0.000 

0.053 

1.094 

1.147 

0.687 

Whitford  Field 

0.865 

0.049  1.343 

2.257 

0.602 

0. 8o5 

0.  052 

1.204 

2.121 

0.  b25 

0.865 

*.053 

1.133 

2.051 

0.625 

NOTE:  i  limited  to  <  a 
e 


Since  variations  in  Ng  in  the  ray  trace  program  have  slight  effect 
on  the  predicted  correlation  bandwidth  and  since  in  the  present  model  the 
scatter  angle  is  based  on  a  4/3  earth  radius,  the  scale  factor  for  the 
condition  of  the  troposphere  (equaling  1  for  turbulent  and  0.4  for  layering) 
is  the  only  factor  which  causes  a  change  in  the  predicted  correlation  band- 
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width  for  a  given  path.  A  method  to  determine  the  condition  of  the  scatter 
volume  from  observations  on  the  surface  was  given  in  Section  II-G  but 
dependable  predictions  will  depend  on  the  availability  of  sonde  data. 

From  these  results  it  can  be  concluded  that  short  term  (1-3  hours) 
variations  in  the  correlation  bandwidth  are  to  be  expected  but  that  they 
cannot  be  accurately  predicted  by  weather  conditions  measured  only  at  the 
surface.  The  predictions  of  the  lower  and  upper  deciles  can  however  be 
made  and  gives  a  statistical  base  for  the  expected  variation  in  the 
correlation  bandwidth. 


p 

e 


Frequency  Separation  —  MHz 


Figure  30.  Change  in  Correlation  Function 
Allowing  Scatter  Angle  to  Become  a  Function 
of  N  .  Point  Petre  Receive  Site. 
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J.  COMPARISON  OF  CORRELATION  BANDWIDTH  MODEL 
WITH  TEST  RESULTS 

To  evaluate  the  correlation  bandwidth  model  presented  in  the  previous 
section,  plots  of  the  correlation  coefficients  predicted  by  the  model  and 
actual  measured  correlation  coefficients  for  eight  different  paths  are 
given  in  Figures  31  through  38.  These  eight  paths  vary  in  length  from 
58  miles  to  200  miles.  Take  off  angles  ranging  from  0.86  degrees  to 
-0.46  degrees  are  included.  Data  from  C-and  X-bands  are  used  for  the 
actual  values. 
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Four  of  the  paths  are  those  used  in  the  test  phase  of  this  program. 

The  effective  beamwidth  function  was  derived  using  the  X-band  lower  deciles 
for  these  four  paths.  There  is  very  good  agreement  between  the  predicted 
correlation  bandwidths  and  the  measured  lower  decile  values  on  these  paths 
(Figures  31  through  34)  .  The  minor  variation  between  the  model  and  experi¬ 
mental  curves  arises  since  the  actual  data  was  not  strictly  Gaussian.  The 
agreement  between  measured  and  predicted  upper  deciles  is  also  good  with 
the  exception  of  Ontario  Center.  This  was  the  shortest  path  with  the 
smallest  scatter  angle  and  suggests  (as  does  the  Tobyhanna-Hexagon  data 
discussed  below)  that  0.4  St  may  fall  to  approximately  0.3  St  for  shorter 
paths.  Alternatively,  it  may  imply  that  the  effective  beamwidth  was  still 
less  than  the  actual  beamwidth  for  this  short  path.  As  discussed  earlier, 
the  values  predicted  by  the  model  are  the  same  for  X-  and  C-band.  These 
predicted  values  can  be  compared  with  C-band  values  plotted  in  the  section 
of  reduced  data. 

C-band  data  taken  over  a  path  from  Verona,  New  York  to  Stony  Point, 

New  York,  in  August  1962  is  plotted  in  Figure  35  along  with  the  values  pre¬ 
dicted  for  this  path^S.  The  lower  decile  predicted  by  the  model  is  about 
one  megahertz  wider  in  correlation  bandwidth  than  the  measured  lower  decile. 
The  measured  upper  decile  is  narrower  in  correlation  bandwidth  than  the 
predicted  upper  decile  using  a  scale  factor  of  0.4 

Figures  36  and  37  have  measured  C-band  data  taken  in  Florida23.  Data 
for  Figure  36  was  obtained  over  the  path  from  Tallahassee  to  Orlando 
during  March  and  April  1963.  The  predicted  lower  decile  correlation  band¬ 
width  agrees  closely  with  the  measured,  but  the  measured  upper  decile 
bandwidth  is  considerably  wider  than  predicted.  Figure  37  has  data 
plotted  which  was  taken  over  a  path  from  Gainesville  to  Orlando  in  May 
and  June  1962.  The  values  predicted  for  this  path  are  much  wider  than 
the  measured  values.  Over  the  Florida  terrain  considerable  difference  in 
the  structure  of  the  scatter  volume  from  the  New  York  area  is  to  be  ex¬ 
pected.  Some  adjustment  of  the  scale  factor  based  on  terrain  could 
improve  the  model  for  paths  of  this  type. 

Tests  run  on  C-band  from  February  through  July  1968  on  the  path  from 
Tobyhanna  to  the  Hexagon  are  plotted  in  Figure  38. 23  The  predicted  values 
are  narrower  than  the  measured  correlation  bandwidths.  It  should  be  pointed 
out,  however,  that  this  path  has  a  small  scatter  angle  and  the  function 
for  effective  beamwidth  in  terms  of  scatter  angle  becomes  inaccurate  or 
at  least  highly  sensitive  for  small  scatter  angles.  Some  refinement  of 
this  function  could  lead  to  better  predictions  for  short  paths  or  for 
paths  with  negative  take  off  angles. 
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Figure  31.  Measured  X-Band,  Summer,  Ontario  Center 
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III.  FIELD  TEST  PROGRAM 

The  purpose  of  the  field  test  program  was  to  record  signal  strength 
data  on  tape  over  several  troposcatter  paths  representative  of  the  MALLARD 
links.  From  the  recorded  tapes  the  cross  correlation  coefficients  versus 
frequency  spacing  and  the  fade  statistics  were  later  determined  through 
computer  processing.  Particular  emphasis  was  placed  on  obtaining  data 
descriptive  of  effects  due  to  terrain,  path  length,  season,  time  of  day, 
metero logical  conditions  and  frequency.  These  results  were  also  used  to 
update  the  analytical  model  for  the  prediction  of  correlation  bandwidth. 
Three  overland  and  one  overwater  path  provided  variations  in  terrain  and 
path  length.  Seasonal  variations  were  obtained  by  changing  paths  every 
few  weeks  over  a  seven  month  testing  period  from  August  1969  through 
February  1970.  The  daily  operating  schedule  was  altered  to  provide  infor¬ 
mation  on  the  diurnal  effects.  Meterological  conditions  were  recorded  and 
are  included  ir.  Appendices  B  and  D.  Throughout  the  program  tests  were 
conducted  simultaneously  on  frequencies  of  4.62  and  7.6  GHz  to  provide 
variations  in  frequency.  Special  tests  were  conducted  over  a  narrow  band¬ 
width  to  determine  the  shape  of  the  correlation  coefficient  for  small 
changes  in  frequency. 

The  transmitters  and  a  ten  foot  diameter  parabolic  antenna  were  located 
at  the  RADC  test  site  in  Model  City,  New  York  (Figure  39).  The  receiving 
instrumentation  was  fully  mobile  and  housed  in  a  32  foot  trailer  with 
special  modifications  to  include  a  ten  foot  stowable  antenna  on  top.  Three 
receiver  sites  were  visisted  cyclically:  Ontario  Center,  New  York,  Whitford 
Field  near  Weedsport,  New  York,  and  Point  Petre ,  Ontario,  Canada.  A  fourth 
site  located  atop  a  hill  at  Port  Byron,  New  York,  was  added  at  the  end  of 
the  program  to  provide  data  comparable  with  the  Whitford  Field  path  but  with 
a  lower  antenna  take  off  angle.  All  the  test  objectives  were  achieved  with 
a  bonus  of  data  from  the  fourth  receiving  site  at  Port  Byron,  New  York. 
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A. 


TRANSMITTER  SITE 


The  transmitter  site  was  located  at  the  RADC  Test  Site  at  Model  City, 

New  York.  The  antenna  (Figure  40)  is  mounted  on  the  roof  of  the  building 
housing  the  transmitter  equipment  and  is  at  an  elevation  of  338  feet  (102 
meters)  above  sea  level.  The  position  of  Model  City  is:  Latitude:  43  deg, 
12  min.,  37.5sec  North;  Longitude:  78  deg,  59  min.,  15.7sec  West. 

The  Ontario  Center  to  Model  City  link  is  a  140  kilometer  path  with  the 
receiver  located  on  a  hill  at  the  RADC  Test  Site  Ontario  Center,  New  York. 
Figure  41  is  a  view  along  the  path  from  Ontario  Center.  Figure  42  is  a 
partial  propagational  profile  chart  showing  the  first  15  miles  from  the 
receive  and  transmit  antennas.  Figure  43  is  a  complete  propagational  profile 
chart  of  the  path.  The  first  obstruction  along  the  path  occurred  5.2  miles 
from  the  receive  antenna  at  an  elevation  of  545  feet  above  sea  level. 


Figure  40.  Antenna 
Installation  at  Model 
City,  N.  Y. 
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Figure  41.  View  Along  Path  From  Ontario 
Center  Toward  Model  City,  N.  Y. 
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Figure  42.  Partial  Profiles  of  Model  City  to  Ontario  Center  Path 
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The  first  obstruction  along  the  transmit  path  occurred  5.8  miles  from  the 
antenna  at  an  elevation  of  380  feet  above  sea  level. 


This  path  is  traversed  by  many  aircraft  entering  and  leaving  the 
Rochester  area.  Aircraft  in  the  path  caused  no  problem  in  reducing  the  data 
as  the  data  reduction  process  was  designed  to  permit  exclusion  of  aircraft 
interference. 


The  weather  conditions  at  Ontario  Center  and  Model  City  were  usually 
similar  to  each  other  throughout  the  test  period. 


Significant  path  parameters  for  this  link  are: 
Position  of  Ontario  Center: 


Latitude : 
Longitude : 
Antenna  Altitude: 


43  degrees,  11  min.,  32sec  North 
77  degrees,  15  min.,  42sec  West 
178.6  meters 


Path  length: 

Bearing  to  Model  City: 

Bearing  to  Ontario  Center: 
Ontario  Center  takeoff  angle: 
Model  City  takeoff  angle: 
Scatter  angle: 

Path  loss: 


86.89  statute  miles,  139.8  km 

271.4  degrees  true 

90.2  degrees  true 

-0.12  degrees 

0.02  degrees 

0.844  degrees 

4.62  GHz,  215.5  dB 
7.6  GHz,  222.9  dB 


The  Point  Petre  to  Model  City  link  is  a  163  kilometer  overwater  path. 

The  receiver  site  was  located  on  the  shore  of  Lake  Ontario  at  Point  Petre' 

Ontario,  Canada  (Figure  44)  .  This  site  afforded  an  unobstructed  view  as 

the  receiver  van  was  parked  only  a  few  feet  from  the  shore.  Figures  45  and 

46  are  partial  and  complete  propagational  profile  charts  of  the  path.  The 

receive  antenna  was  at  an  elevation  of  245  feet  above  sea  level  with  the 

path  grazing  the  water  approximately  5  miles  from  the  antenna.  The  first 

obstruction  along  the  transmit  path  occurred  7.6  miles  from  the  antenna  at 

an  elevation  of  340  feet  above  sea  level. 

i 

r 

Aircraft  interference  was  rarely  observed  on  this  path. 

The  weather  conditions  at  Point  Petre  and  Model  City  were  similar  to 
each  other  during  the  summer  testing.  However,  during  the  winter  (December 
and  January)  there  were  days  with  noticeable  differences  in  weather  condi¬ 
tions.  The  general  weather  condition  for  these  days  was  a  northerly  wind 
with  the  temperature  at  Point  Petre  10  to  25  degrees  colder  than  Model  City. 
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Figure  44.  Point  Petre  Receiving  Site 
Significant  parameters  for  this  link  are: 


Position  of  Point  Petre  : 

Latitude : 
Longitude : 

Antenna  Altitude: 

Path  length : 

Bearing  to  Model  City: 

Bearing  to  Point  Petre  : 

Point  Petre  takeoff  angle: 

Model  City  takeoff  angle: 

Scatter  angle: 


43  degrees,  50  min.,  22sec  North 
77  degrees,  9  min.,  16sec  West 
80.77  meters 

101.5  statute  miles,  163.4  km 

245.3  degrees  true 

64.0  degrees  true 

-0.06  degrees 

-0.03  degrees 

1.011  degrees 


Path  loss:  4.62  GHz,  219.5  dB 

7.6  GHz,  227.1  dB 


The  Whitford  Field  to  Model  City  link  is  a  199  kilometer  overland  path. 
The  receive  site  was  located  at  Whitford  Field  near  Weedsport,  New  York. 
Figure  47  is  an  aerial  photograph  of  the  site  with  a  major  obstruction  along 
the  path  marked.  This  was  a  marginal  path  as  there  were  two  major  obstruc¬ 
tions  near  the  antenna.  The  first  obstruction  was  a  hedgerow  1,500  feet 
from  the  antenna.  This  obstruction  was  eliminated  by  cutting  a  path  through 
the  trees  after  the  antennas  were  aligned  in  azimuth.  The  second  obstruc- 
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Figure  46.  Metric  Profile  of  Model  City  to  Point  Petre  Path 


Figure  47.  Whitford  Field 


tion  was  also  caused  by  trees  2,000  feet  from  the  antenna.  A  path  could  not 
conveniently  be  cut  through  these  trees  as  there  were  too  many.  The 
trees  caused  the  takeoff  angle  to  be  0.862  degrees  or  0.442  degrees  greater 
than  first  calculated.  Figures  48  and  49  are  partial  and  complete  propaga- 
tional  profile  charts  oi  the  path.  The  receive  antenna  was  at  an  elevation 
of  410  feet  above  sea  level  with  the  first  obstruction  occurring  2,000  feet 
from  the  antenna  at  an  elevation  of  440  feet  above  sea  level.  The  first 
obstruction  along  the  transmit  path  occurred  5.0  miles  from  the  antenna  at 
an  elevation  of  370  feet  above  sea  level. 

Aircraft  interference  on  this  path  was  present  5  to  10  percent  of  the 
time.  The  center  of  the  beam  passed  within  one  mile  of  the  Rochester 
airport  and  the  lower  ray  was  at  an  altitude  of  2,700  feet  above  sea  level. 

The  weather  conditions  at  Whitford  Field  and  Model  City  were  usually 
similar  throughout  the  test  period. 

Significant  path  parameters  for  this  link  are: 


51 


PROFILE 

BASED  ON  1.33  x  EARTH'S  RADIUS 


m  rfpniisiv 


52 


Lgure  48.  Partial  Profile  of  Model  City  to  Whitford  Field  Path 
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Figure  49.  Metric  Profile  of  Model  City  to  Whitford  Field  Path 


Position  of  Whitford  Field: 

Latitude:  43  degrees,  4  min.,  45sec  North 

Longitude;  76  degrees,  32  min.,  22sec  West 

Antenna  Mtitude: 

Path  length:  123.7  statute  miles,  199.2  km 

Bearing  to  Model  City:  275.0  degrees  true 

Bearing  to  Whitford  Field:  93.3  degrees 

Whitford  Field  takeoff  angle:  0.862  degrees 

Model  City  takeoff  angle:  0.05  degrees 

Scatter  angle:  2.15  degrees 

Path  loss:  4.62  GHz,  232.1  dB 

7.6  GHz,  240.0  dB 

The  Port  Byron  to  Model  City  link  is  a  192  kilometer  path  with  the 
receiver  van  located  on  top  of  a  hill  in  Port  Byron,  New  York  (Figure  50)  . 

This  site  was  used  to  obtain  comparative  data  over  a  path  similar  to  the 
Whitford  Field  to  Model  City  path.  Figures  51  and  52  are  partial  and  complete 
propagational  profile  charts  of  the  path.  The  receive  antenna  was  at  an 
altitude  of  564  feet  above  sea  level  with  the  first  obstruction  along  the 
receive  path  occurring  12  miles  from  the  antenna  at  an  elevation  of  620  feet 
above  sea  level.  The  first  obstruction  along  the  transmit  path  occurred 
6.2  miles  from  the  antenna  at  an  elevation  of  392  feet  above  sea  level. 

Aircraft  interference  was  more  numerous  on  this  path  than  on  any  of  the 
other  paths.  Aircraft  interference  though  eliminated  in  data  reduction  was 
present  in  30  percent  of  the  raw  data.  The  center  of  the  beam  passed 
directly  over  the  Rochester  airport,  and  the  lower  ray  was  at  an  altitude 
of  1900  feet  above  sea  level.  The  intersection  of  the  lower  rays  occurred 
8.1  miles  west  of  the  airport. 

Complete  weather  information  at  Port  Byron  was  not  available  as  the 
closest  weather  services  were  located  at  Syracuse  (40  miles)  and  Rochester 
(50  miles) .  The  Rochester  weather  service  was  used  since  the  common  volume 
was  close  to  the  weather  bureau. 

Significant  path  parameters  are  listed  below: 

Position  of  Port  Byron: 

Latitude:  43  degrees,  2  min.,  2.4sec  North 

Longitude:  76  degrees,  37  min.,  51.6sec  West 

Antenna  Altitude:  184  meters,  564  feet 
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Position  of  Whitford  Field: 


Latitude:  43  degrees,  4  min.,  45sec  North 

Longitude:  76  degrees,  32  min.,  22sec  West 

Antenna  Altitude: 

Path  length:  123.7  statute  miles,  199.2  km 

Bearing  to  Model  City:  275.0  degrees  true 

Bearing  to  Whitford  Field:  93.3  degrees 

Whitford  Field  takeoff  angle:  0.862  degrees 

Model  City  takeoff  angle:  0.05  degrees 

Scatter  angle:  2.15  degrees 

Path  loss:  4.62  GHz,  232.1  dB 

7.6  GHz,  240.0  dB 

The  Port  Byron  to  Model  City  link  is  a  192  kilometer  path  with  the 
receiver  van  located  on  top  of  a  hill  in  Port  Byron,  New  York  (Figure  50) . 

This  site  was  used  to  obtain  comparative  data  over  a  path  similar  to  the 
Whitford  Field  to  Model  City  path.  Figures  51  and  52  are  partial  and  complete 
propagational  profile  charts  of  the  path.  The  receive  antenna  was  at  an 
altitude  of  564  feet  above  sea  level  with  the  first  obstruction  along  the 
receive  path  occurring  12  miles  from  the  antenna  at  an  elevation  of  620  feet 
above  sea  level.  The  first  obstruction  along  the  transmit  path  occurred 
6.2  miles  from  the  antenna  at  an  elevation  of  392  feet  above  sea  level. 

Aircraft  interference  was  more  numerous  on  this  path  than  on  any  of  the 
other  paths.  Aircraft  interference  though  eliminated  in  data  reduction  was 
present  in  30  percent  of  the  raw  data.  The  center  of  the  beam  passed 
directly  over  the  Rochester  airport,  and  the  lower  ray  was  at  an  altitude 
of  1900  feet  above  sea  level.  The  intersection  of  the  lower  rays  occurred 
8.1  miles  west  of  the  airport. 

Complete  weather  information  at  Port  Byron  was  not  available  as  the 
closest  weather  services  were  located  at  Syracuse  (40  miles)  and  Rochester 
(50  miles) .  The  Rochester  weather  service  was  used  since  the  common  volume 
was  close  to  the  weather  bureau. 

Significant  path  parameters  are  listed  below: 

Position  of  Port  Byron: 

Latitude:  43  degrees,  2  min.,  2.4sec  North 

Longitude:  76  degrees,  37  min.,  51.6sec  West 

Antenna  Altitude:  184  meters,  564  feet 
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Figure  50.  Port  Byron  Site 


Path  Length : 

Bearing  to  Model  City: 
Bearing  to  Port  Byron: 
Model  City  takeoff  angle: 
Port  Byron  takeoff  angle: 


Scatter  angle: 


Path  loss: 


119  statute  miles,  192  km 


266.7  degrees  true 


95.0  degrees  true 


0.49  degrees 

-0.083  degrees 

1.335  degrees 

4.62  GHz,  222.4  dB 
7.6  GHz,  230.2  dB 
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Lgure  51 .  Partial  Profile  of  Model  City  to  Port  Byron  Path 


Figure  52.  Metric  Profile  of  Model  City  to  Port  Byron  Path 


B.  TRANSMITTER  INSTRUMENTATION 

The  transmitting  instrumentation  was  located  in  the  RADC  building  at 
Model  City,  New  York.  X-band  signals  at  7.6  GHz  and  1000  watts  were  provided 
by  an  FRC-68  transmitter  which  was  permanently  installed  in  the  building. 

The  C-band  signals  at  500  watts  were  provided  by  an  AN/GRC-66  transmitter 
with  power  amplifier.  Two  types  of  modulation  were  used  to  provide  output 
signals  that  were  effectively  multiple  carrier.  A  system  of  cyclically 
pulsing  two  sets  of  five  gated  oscillators  with  a  100  microsecond  dwell 
on  each  frequency  was  used  to  yield  signals  that  were  incrementally 
separated  from  1  to  9  MHz.  An  FM  signal  was  also  used  at  a  modulation 
index  of  1.841  and  200  kHz  to  enable  correlation  coefficient  measurements 
for  small  increments  in  frequency. 

The  transmitting  equipment  (Figure  53)  transmits  on  both  C-  and  X-bands 
simultaneously.  In  the  wideband  mode,  the  oscillator  stepper  gates  each  of 
the  five  oscillators  cyclically  into  the  70  MHz  input  of  each  transmitter. 

The  oscillator  frequencies  for  the  X-band  modulation  pattern  are  65.5,  66.5, 
69.5,  72.5,  and  74.5  MHz.  This  separation  after  amplification  and  transla¬ 
tion  in  the  FRC-68  transmitter  results  in  stepped  carriers  separated  by  1, 

4,  7,  and  9  MHz  from  the  lower  frequency.  The  stepped  carriers  have  a  peak 
power  of  1000  watts  and  a  duty  cycle  of  0.2.  This  method  of  pulsing  results 
in  only  one  frequency  being  amplified  by  the  transmitter  power  amplifier  at 
a  time  and  thus  there  are  no  intermodulation  products  generated. 
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Figure  53.  Transmitting  Equipment 
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The  C-band  AN/GRC-66  transmitter  goes  through  one  doubling  process  prior  to 
the  final  output.  Thus  a  second  set  of  oscillators  having  separations  half 
that  of  the  other  equipment  was  provided.  These  frequencies  were  67.75, 
68.25,  69.75,  71.25,  and  72.25  MHz.  A  spectrum  analyzer  was  used  in  each 
system  for  monitoring  the  RF  output  spectrum  and  as  an  aid  in  adjusting  the 
equipments.  The  output  signals  of  the  two  transmitter  power  amplifiers  were 
connected  through  waveguide  to  a  dual  frequency  feed  to  illuminate  a  10  foot 
parabola  mounted  on  the  roof  of  the  building  (Figure  54) . 

The  narrow  band  mode  used  the  standard  FM  inputs  of  both  transmitters. 
These  were  set  to  a  modulating  frequency  of  200  kHz  at  an  index  of  1.841  to 
produce  the  set  of  five  significant  spectral  lines  separated  by  200  kHz  and 
having  amplitude  differences  of  about  5  dB. 
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Figure  54.  Antenna  with  Dual  Frequency  Feed 
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The  antenna  characteristics  are  given  in  the  table  below: 


Table  III.  ANTENNA  CHARACTERISTICS 


Antenna 

Diameter  -  10  feet 


Frequency 

3  dB  BW 

(deg) 

Gain 

First 

Side  Lobe 

(dB) 

(GHz) 

E 

H 

(dB) 

VSWR 

E 

H 

4.62 

1.5 

1.4 

41.2 

1.28 

-26 

-20 

7.60 

1.07 

0.92 

43.9 

1.04 

-26 

-19 

Note:  Tests  conducted  at  the  Martin  Marietta  Antenna  and  Microwave 

Laboratory  showed  no  deviations  from  concentricity  between 
the  C-  and  X-band  beams . 


C.  RECEIVING  INSTRUMENTATION 

All  receiver  and  data  recording  equipment  necessary  to  operate  the 
receiver  site  was  mounted  in  a  32-foot  van  (Figure  55) .  The  receiver  site 
block  diagram  is  shown  in  Figure  56.  A  10-foot  parabolic  reflector  with  a 
dual  frequency  feed,  similar  to  a  dual  frequency  transmitter  feed  was 
mounted  on  the  roof  of  the  van  with  its  structural  support  through  the 
trailer  to  the  frame.  The  dish  was  adjustable  between  -3  and  +7  degrees  in 
elevation  and  360  degrees  in  azimuth.  In  transit,  the  feed  was  removed  and 
the  antenna  was  lowered  in  a  horizontal  position.  The  dual  frequency  feed 
permitted  simultaneous  reception  at  4.62  and  7.6  GHz.  The  AN/GRC-66  and  the 
FRC-68  receivers  were  adjusted  to  an  overall  bandwidth  of  10  MHz.  No  AGC 
action  takes  place  in  the  units. 

For  signal  strength  measurements  with  the  wideband  frequency  spacings, 
the  70  MHz  IF  outputs  of  the  receivers  were  connected  through  two  5-port 
multicouplers  to  10  individual  IF  amplifiers.  These  amplifiers  were  tuned 
to  center  frequencies  of  65.5,  66.5,  69.5,  72.5  and  74.5  MHz.  Signal  strength 
measurements  using  the  narrow  frequency  spacings  required  a  second  conversion 
of  the  70  MHz  output.  This  was  accomplished  by  mixing  the  70  MHz  output  with 
a  63.5  MHz  crystal  oscillator  to  give  an  output  around  6.5  MHz.  These 
signals  were  routed  through  a  5-port  multicoupler  to  five  R390A-URR  general 
purpose  receivers  with  bandwidths  of  10  kHz  and  center  frequencies  which 
were  easily  adjusted  to  6.10,  6.30,  6.50,  6.70,  and  6.90  MHz. 

The  AGC  voltages,  which  are  indications  of  signal  strength,  from  the 
10  IF  amplifiers  and/or  the  five  separate  R390  receivers  were  routed  through 
15  operational  amplifiers.  The  amplifiers  were  adjusted  so  that  the  outputs 
varied  between  ±2.00  volts  for  a  50  dB  change  in  signal  strength. 
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Figure  55.  Receiving  Instrumentation  Van 


Ten  of  the  signal  conditioner  outputs  are  routed  through  the  automatic 
test  identity  and  control  panel  to  60  Hz  notch  filters.  These  filters  remove 
any  60  Hz  residual  hum  in  the  system.  These  active  notch  filters  are  only  a 
few  Hertz  wide  at  60  Hz.  Outputs  from  the  notch  filters  are  inserced  into  the 
FM  recording  oscillators  which  are  part  of  the  0.6  ips  tape  recorder.  These 
oscillators  are  linear  in  their  input/output  characteristics.  The  FM  record 
oscillators  v/ere  adjusted  so  that  an  input  between  ±2.00  volts  caused  the 
output  to  vary  between  740  and  340  Hz.  Parallel  outputs  from  the  notch 
filters  also  go  to  a  strip  chart  recorder  that  provides  a  visual,  on-site 
estimation  of  fade  rates  and  a  check  on  system  performance. 

The  automatic  test  identity  and  control  panel  was  provided  to  insert 
the  test  number  and  weather  data  directly  on  the  tape  to  avoid  any  discrep¬ 
ancies  between  test  logs,  etc.  The  control  panel  also  automatically 
controlled  the  sequence  of  events  from  start  to  finish  of  each  test. 

The  overall  receiver  to  tape  recorder  calibration  curve  was  provided  to 
the  Martin  Marietta  computer  so  that  the  resulting  data  could  be  related  to 
real  time  signal  strengths  in  dBm.  The  narrow  band  calibration  was  performed 
with  a  70  MHz  signal  generator  to  generate  CW  signals  of  known  signal 
strength  into  the  70  MHz  input  of  the  narrow  band  system.  Calibration  of  the 
wideband  receiving  system  was  accomplished  using  a  pulsed  signal  generator 
at  4.62  or  7.6  GHz  at  the  RF  input  to  the  receivers. 


IF  Amplifiers 
\\  ideband 
Frequencies 

t  AGC  Lines 


Operational  Amplify 
I 


"X"  Band 
Waveguide 


Narrow  Band 
Receivers 
(R390's) 


"C"  Band 
Waveguide 


74.5  M  Hz 

F5 

"1 

IF 

_ 

72.5  MHz 

F  4 

IF 

_ 

69.5  MHz 

F3 

IF 

_ 1 

66.5  MHz 

F2 

- 1 

IF 

_ \ 

65.5  MHz 

FI 

1 

IF 

_ r 

H 


F  5  __ 

6.7  MHz 

F4  _ 

F3  ^ 

6.3  MHz 

F2« 

F1  • 

u ,  i  ivi  n  z 

74.5  MHz 

F5 

1 - 1 

t 

1  1 

72.5  MHz 

I  F  4 

* 

1 

IF 

_ 

69.5  MHz 

F3 

IF 

66.5  MHz 

F2 

1 

IF 

1 

65.5  MHz 

FI  . 

IF 

H 


i 


Figure  56.  Receiving  Instrumentation 
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D.  ANTENNA  ALIGNMENT 

The  receive  antenna  was  aligned  to  the  transmit  antenna  by  a  coarse 
boresight  procedure  and  by  the  use  of  the  maximum  of  maxima  method.  The 
first  step  in  the  alignment  was  to  use  a  hand  compass  at  both  receiver  and 
transmitter  sites  to  obtain  coarse  azimuth  settings.  The  boresight  was  then 
used  for  the  coarse  elevation  settings.  The  accuracy  of  these  settings  was 
usually  +3  degrees  in  azimuth  and  +1.5  degrees  in  elevation.  The  transmitter 
was  then  set  up  to  transmit  a  7.6  GHz  CW  signal  with  a  69.5  MHz  output  from 
the  receiver.  The  receive  antenna  was  then  swept  in  the  azimuth  plane  across 
the  transmitted  beam.  A  strip  chart  recording  of  the  69.5  MHz  IF  amplifier 
AGC  was  made  during  the  sweep.  The  receive  antenna  was  then  set  at  the 
position  where  the  strongest  received  signal  occurred  and  the  position  of  the 
antenna  was  logged.  The  same  procedure  was  then  performed  on  the  transmit 
antenna.  The  overall  procedure  served  to  align  the  antennas  within  +0.75  degrees 
of  final  alignment.  The  average  signal  strength  was  measured.  The  receive  antenna 
was  then  swept  across  the  beam  in  steps  of  0.3  degrees.  The  antenna  position  and 
average  signal  strength  was  recorded  for  each  step.  The  transmit  antenna  was 
then  moved  0.5  degrees  to  the  left  or  the  right.  The  receive  antenna  is  again 
swept  through  the  beam  in  0.3  degree  steps  with  antenna  position  and  average 
signal  strength  recorded.  This  process  was  repeated  in  smaller  increments 
until  there  was  no  further  improvement  in  received  signal.  The  transmit  and 
receive  antennas  were  then  locked  down  in  the  azimuth  plane.  The  same  proce¬ 
dure  was  used  for  the  elevation  plane. 
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IV.  DATA  COLLECTION 

During  the  testing  period  from  August  1969  through  February  1970  data 
were  collected  over  the  troposcatter  paths.  The  signal  strength  data  were 
recorded  on  magnetic  tape  covering  both  C-  and  X-bands  simultaneously  during 
the  conventional  wideband  tests.  The  special  narrow  bandwidth  tests  were 
made  singly  since  only  one  set  of  five  narrowband  R390  communications 
receivers  were  available  for  the  program  and  10  receivers  would  have  been 
required  for  simultaneous  measurements.  Supporting  data  were  also  required. 
Weather  records  were  obtained  from  the  Weather  Bureau  at  Rochester,  New  York. 
These  are  included  in  Appendix  B.  Weather  data  during  each  test  was  also 
obtained  on  a  real-time  basis  by  telephone  from  the  Buffalo,  New  York, 
Rochester,  New  York,  and  the  Trenton,  Ontario  Weather  Bureaus.  The  real 
time  weather  data  was  inserted  on  the  tapes  and/or  included  in  the  operator's 
test  log.  These  data  are  included  among  other  data  in  Appendix  D,  Test 
Records.  Average  signal  strength  was  obtained  from  a  counter  whose  input 
from  one  of  the  FM  record  oscillators  represented  the  time  varying  signal 
strength.  This  figure  was  entered  in  the  operator's  log.  The  actual  median 
signal  strength  was  obtained  through  the  computer  reduction.  The  median 
value  of  each  test  reduced  is  posted  in  Appendix  D. 

The  schedule  of  operation  at  the  various  sites  was  planned  to  provide  a 
representation  of  seasonal  variations  as  well  as  terrain  variations.  Table 
IV  is  a  list  of  the  sites  and  the  testing  periods  at  each.  The  last  site 
listed  is  Port  Byron  which  is  a  site  similar  in  path  length  to  the  Whitford 
Field  site,  with  the  exception  of  a  lower  antenna  takeoff  angle. 


Table  IV.  SITE  OPERATION  SCHEDULE 


Ontario  Center ,  Summer 
Whitford  Field,  Summer 
Point  Petre,  September 
Ontario  Center,  October 
Whitford  Field,  November 
Point  Petre,  Winter 
Ontario  Center,  Winter 
Port  Byron,  February 


August  4,  1969,  through  August  20,  1969 
August  26,  1969,  through  September  9,  1969 
September  15,  1969,  through  September  26,  1969 
October  1,  1969,  through  October  24,  1969 
October  29,  1969,  through  November  18,  1969 
December  2,  1969,  through  January  15,  1970 
January  22,  1970,  through  February  5,  1970 
February  7,  1970,  through  February  26,  1970 


The  daily  operating  schedule  was  rotated  in  order  to  assess  diurnal 
effects . 
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The  receiving  systems  were  calibrated  as  often  as  necessary  to  ensure 
that  the  signal  strength  data  were  accurate  to  the  tolerance  of  the  measuring 
equipment.  Accuracies  of  about  one  dB  are  typical  except  at  signal  strengths 
near  the  receiver  thresholds.  The  signal  strength  itself  was  not  used  as  a 
variable.  Its  median  value,  however,  was  logged  (Appendix  D)  and  used  to 
determine  the  sufficiency  of  fade  margin  to  ensure  the  validity  of  the  depth 
of  fade  and  signal  amplitude  distributions. 

Except  as  noted  in  the  data  and  test  logs  the  power  output  of  the  X-band 
transmitter  was  1000  watts  peaK  with  a  duty  cycle  of  0.2  on  the  wide 'band 
tests  and  the  C-band  tests  used  500  watts  and  the  same  duty  cycle.  The 
special  narrow  bandwidth  tests  were  using  FM  at  1000  watts  for  X-band  and 
500  watts  for  C-band. 

Once  each  week  the  magnetic  tapes  were  mailed  from  the  operating  site  to 
the  Orlando  facility  for  processing  through  the  data  reduction  facility.  The 
actual  data  processing  is  described  in  the  following  section  of  this  volume. 


V. 


DATA  PROCESSING 


A.  INTERMEDIATE  PROCESSING 

Signal  strength  data  for  the  correlation  bandwidth  study  was  recorded 
on  ten  channels  of  an  analog  tape  at  0.6  ips  as  previously  described.  A 
strip  chart  recording  of  either  X-  or  C-band  data  was  also  made.  A  list  of 
tests  with  corresponding  information  about  bandwidth,  AGC  calibrations, 
weather,  and  special  test  conditions  was  kept  in  the  test  log. 

When  the  recorded  tapes  were  received  at  the  Martin  Marietta  plant, 
they  were  processed  by  MADRE,  Martin  Automatic  Data  Reduction  Equipment. 

MADRE  processed  the  tape  separately  for  X-  and  C-band  signals.  Each  time 
five  channels  were  digitized  and  recorded  along  with  a  time  channel  on  an 
intermediate  tape  (see  Figure  57)  .  A  Brush  strip  chart  recording  of  the 
five  data  channels  and  a  time  marker  corresponding  to  the  time  channel  on 
the  intermediate  tape  also  was  made.  C-band  data  recorded  on  channels  1,  3, 
5,  7,  and  9  was  played  back  at  a  50:1  speed  up.  X-band  data  from  channels 
2,  4,  6,  8,  and  12  was  digitized  at  a  speed  up  of  25:1.  The  lower  speed  up 
was  necessary  to  ensure  adequate  sampling  of  the  higher  fade  rates  at  X-band. 
MADRE' s  sample  rate  of  20,000  samples  per  second  resulted  in  a  sampling- of 
each  X-band  channel  every  1/160  second  of  actual  test  time  and  each  C-band 
channel  every  1/80  second  of  real  test  time. 

After  the  digitizing  was  completed,  the  Brush  charts,  strip  charts 
made  at  the  receiver  site,  and  test  logs  were  used  to  select  parts  of  the 
tests  where  all  channels  were  operating  correctly  and  which  were  not  affected 
by  aircraft.  Some  samples  were  chosen  specifically  to  include  aircraft 
fading  so  its  effect  could  be  studied  separately.  From  acceptable  tests,  a 
segment  of  6  MADRE-time  seconds  was  selected  to  be  reduced  by  the  computer 
program.  For  C-band  with  a  50:1  speed-up  6  MADRE  seconds  equals  five  minutes 
of  field  testing.  On  X-band  with  a  25:1  speed-up  the  6  MADRE  seconds  repre¬ 
sented  two  and  a  half  minutes  of  field  testing.  The  selected  parts  of  tests 
were  reformatted  and  recorded  on  a  tape  suitable  for  CDC  6400  computer  input 
(see  Figure  58). 

Information  obtained  from  the  test  log  was  prepared  for  the  computer 
data  reduction.  A  preliminary  program  was  used  on  the  CDC  6400  digital 
computer  to  process  the  AGC  calibrations  (Figure  59)  .  Punched  cards  describ¬ 
ing  the  AGC  calibration  curves  were  the  output  from  this  program.  For  each 
test  an  identity  card  was  manually  punched  for  both  X-  and  C-band. 
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Figure  57.  Analog  to  Digital  Conversion 
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Figure  58.  Final  Digital  Formatting  of  Selected  Samples 
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This  card  contained  information  about  the  test  including  date,  time,  test 
number,  frequency,  bandwidth,  receive  site  weather  and  temperature,  trans¬ 
mit  site  weather  and  temperature,  and  AGC  curve  identification.  The  test 
identification  word  was  used  on  all  the  printed  and  punched  data  for  the 
test.  A  typical  test  identification  can  be  interpreted  as  follows: 

02  17  1320  09  W  C  21  16  05  16 

Transmit  Site  Temperature  Code 
Transmit  Site  Weather  Code 
Receive  Site  Temperature  Code 
Receive  Site  Weather  Code 
C-Band 

Wide  Frequency  Spacing 
Ninth  Test  of  the  Day 
Test  Started  at  1320 
Day 
Month 

The  test  identification  cards  and  AGC  curve  cards  were  input  data  to 
the  computer  program  used  to  read  the  tapes  prepared  by  MADRE. 

B.  DATA  REDUCTION 

After  MADRE  completed  reformatting  the  data  onto  tapes  for  the  CDC 
6400,  the  tapes  were  read  by  a  program  entitled  SPEED.  The  input-output 
and  some  arithmetic  statements  of  the  program  are  written  in  FORTRAN, 
but  most  of  the  calculations  are  controlled  by  subroutines  written  in 
COMPASS,  a  CDC  machine  language. 

A  separate  run  of  SPEED  was  made  for  each  reformatted  tape  as  in 
Figure  60.  In  addition  to  the  tape  input,  test  identifications  and  AGC 
curves  corresponding  to  the  tests  on  the  tape  were  input  in  card  form. 

The  program  read  the  tape  in  blocks  of  2000  samples  and  calculated  the 
correlation  coefficients,  median  signal  strength,  and  mean  signal 
strength.  Nine  correlation  coefficients  were  calculated  for  wide 
spaced  tests  since  nine  different  frequency  spacings  were  obtained 
by  selected  frequency  separations. 
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Printed  and  Punched  Data 
from  Tape  1 


Program  Speed,  AGC  Curves  and 
Test  Identity  Cards  for  Tests  on 
Tape  1 


Program  Speed,  AGC  Curves  and 
Test  Identity  Cards  for  Tests  on 
Tape  2 


Printed  and  Punched  Data 
from  Tape  2 


Figure  60.  Tape  and  Card  Input  Data  Processing  by  Main  Data 

Reduction  Program 
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Wide  FI  F2 _ F3 _ F4 _ F5 

Spacing  0123456789  MHz 

p  '  nurow  spaced  tests,  only  four  correlation  coefficients  were  calcu- 
•  >  :  mce  the  five  channels  were  uniformly  spaced.  Further  combinations 

ul  t.r.e  channels  would  have  given  only  duplicate  frequency  separations  be¬ 
tween  the  correlated  channels. 

Narrow  Fl _ F2 _ F3 _ F4 _ F5 

Spacing  0  100  200  300  400  500  600  700  800  kHz 

The  decision  to  calculate  nine  or  four  correlation  coefficients  was  made 
by  the  program  based  on  the  letter  W  for  wide  or  N  for  narrow  in  the  test 
identity. 

Fade  rate,  fade  depth,  fade  duration,  and  signal  amplitude  dis¬ 
tribution  calculations  were  made  for  only  one  channel,  F3.  For  each  test, 
all  these  data  were  outputted  by  both  the  printer  and  the  card  puncher.  In 
addition  to  calculated  data,  each  punched  card  contained  the  test  identi¬ 
fication  for  the  test  and  a  code  punch  in  columns  79  and  80  to  identify 
the  type  of  data  on  the  card.  The  printed  outputs  from  the  program  runs 
were  used  to  determine  whether  the  data  was  processed  correctly  and  to 
check  characteristics  of  individual  tests. 

C.  PLOTTING 

Plots  of  the  data  for  some  tests  were  made  either  separately  or 
with  a  number  of  tests  grouped  on  the  same  axis.  All  plot  programs  were 
first  run  on  the  CDC  6400  from  which  a  magnetic  tape  was  made  for  use  on 
the  CALCOMP  763  off-line  plotter  and  CALCOMP  770  console  combination  (see 
Figure  61) . 

When  all  the  data  tapes  for  a  testing  period  at  a  site  had  been 
read  by  SPEED,  the  data  was  machine  sorted  by  the  codes  representing 
frequency,  date  and  type  of  data  desired  (Figure  62).  Prior  to  final 
collation  of  the  data  tor  a  testing  period,  the  printouts  for  all  the 
tests  were  manually  checked.  Tests  for  which  the  median  signal  strengths 
were  not  at  least  20  dB  above  noise  were  considered  to  have  insufficient  fade 
margin  to  be  considered  in  the  fade  depth,  fade  duration,  and  signal  amplitude 
summaries.  Fade  rate,  however,  considers  only  positive  crossing  of  the  median, 
thus  insufficient  fade  margin  had  no  adverse  effects  on  the  fade  rates 
observed.  A  list  of  all  tests  from  which  samples  were  selected  was  pre¬ 
pared  by  the  computer.  This  test  list  included  the  test  number,  receiver 
location,  weather  conditions  at  both  sites,  and  the  median  signal  strength 
during  the  test. 

Distribution  plots  were  made  to  display  the  data  obtained  during 
each  period  of  testing  at  the  various  sites.  Five  kinds  of  data  were 
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plotted  for  both  X-and  C-band.  These  were  correlation  coefficient,  fade 
rate,  fade  depth,  fade  duration,  and  signal  amplitude.  Diurnal  effects 
were  considered  by  making  distribution  plots  for  the  data  of  tests  started 
during  four  time  classes:  0001  to  0600,  0601  to  1200,  1201  to  1800,  and 
1801  to  2400  hours.  Separate  distribution  plots  were  also  made  for  the 
data  of  all  tests  run  when  the  receiver  site  temperature  was  greater  than 
an  estimated  average  temperature  for  the  period  of  testing  and  when  the 
temperature  was  less  than  or  equal  to  the  average  temperature.  In  addition, 
a  distribution  plot  designated  ALL  was  made  for  the  data  of  all  tests  per¬ 
formed  during  the  test  period  without  regard  to  time  of  day  or  tenperature. 

Sets  of  percentile  distributions  were  made  for  correlation  coef¬ 
ficient,  fade  rate,  fade  depth,  and  signal  amplitude  data.  If  ten  or 
more  tests  were  included  in  the  distribution,  the  upper  decile,  median, 
and  lower  decile  were  plotted.  If  fewer  than  three  tests  were  available 
for  a  distribution,  no  plot  was  made.  For  distributions  including  three 
to  nine  tests,  only  the  median  was  plotted. 

In  order  to  summarize  fade  duration  data,  a  cumulative  array  of 
fade  durations  for  each  of  five  levels  from  0  to  20  dB  below  the  median 
was  made.  The  cumulative  percentage  by  fade  duration  class  was  then  cal¬ 
culated  at  each  level  and  plotted. 

Summary  bar  graphs  were  made  to  permit  an  easy  comparison  of  the 
data  obtained  during  different  periods  of  testing  at  a  site.  Separate 
graphs  were  made  to  X-and  C-band  data  from  all  four  sites.  The  infor¬ 
mation  contained  in  the  summary  bar  graph  was  read  from  decile  distribution 
plots.  For  the  correlation  coefficient  data,  the  value  of  the  frequency 
separation  of  the  median  and  upper  and  lower  deciles  for  which  the  cor¬ 
relation  coefficient,  p  ,  equals  0.4  on  the  distribution  plot  was 
plotted  in  the  corresponding  location  on  the  bar  chart.  For  the  fade  rate 
and  fade  depth  data,  the  value  in  Hertz  or  dBm  below  the  median  where  the 
median  and  upper  and  lower  decile  crossed  50  percent,  was  transferred  to 
the  bar  graph.  In  the  signal  amplitude  distributions,  the  value  exceeded 
95  percent  of  the  time  was  plotted  on  the  summary  bar  graph.  In  cases 
where  fewer  than  ten  tests  were  included  in  the  distribution,  only  the 
median  was  read  from  the  plot.  The  fade  duration  summary  bar  graph  was 
done  in  a  similar  manner  but  values  were  read  for  0,  10,  and  20  dB  below  the 
median. 
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VI .  REDUCED  DATA 


In  this  section  distribution  and  summary  plots  of  data  taken  during 
eight  different  test  periods  at  the  four  sites  are  shown.  Correlation 
coefficient,  fade  rate,  fade  depth,  fade  duration,  and' signal  amplitude 
plots  are  included. 

A.  FREQUENCY  CROSS-CORRELATION  MEASUREMENTS 
1.  Individual  Tests  and  Narrow  Spacing  Tests 

The  results  of  individual  tests  made  simultaneously  at  X-  and  C-bands 
are  of  significant  interest.  The  study  indicates  that:  1)  Correlation  band¬ 
width  is  a  rather  rapidly  changing  function  of  time.  2)  There  is  no 
significant  loss  of  frequency  diversity  expected  as  a  result  of  employing 
X-band  as  compared  to  C-band  over  these  paths  and  with  these  antenna 
sizes.  3)  The  correlation  bandwidth  shows  some  relation  with  fade  rate, 
but  not  always  the  same  trend  is  displayed.  4)  For  the  narrow  spacing 
tests  the  derivative  of  the  cross-correlation  coefficient  versus  fre¬ 
quency  separation  is  zero  as  the  frequency  separation  approaches  zero. 

The  data  from  individual  tests  are  plotted  in  Figures  63  through  79. 

A  number  of  individual  tests  are  plotted  on  the  same  axis  to  permit  a 
comparison  of  results  from  one  test  to  the  next.  The  test  numbers  in¬ 
dicated  on  the  curves  are  the  individual  test  identification  numbers 
which  can  be  interpreted  as  follows: 

09  25  1118  06 

TTTL  Sixth  Test  of  the  Day 

Test  Started  at  1118 
The  TVenty  Fifth  Day  of  the  Month 
■  The  Ninth  Month 

Figures  63  and  64  represent  a  set  of  simultaneous  measurements  of 
correlation  versus  frequency  separation  taken  at  C-  andX-bands.  It  is 
immediately  apparent  that  we  are  dealing  with  a  rapidly  changing  phen- 
nomencn  since  from  one  test  to  another  no  two  curves  coincide.  Moreover, 
the  marked  break  between  1050  and  1115  hours  suggests  that  the  scattering 
mechanism  went  through  a  complete  and  abrupt  change  in  its  characteristics. 
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In  Figure  64,  the  same  type  of  phenomena  is  present  at  X-band  and  appar¬ 
ently  occurred  during  the  same  time  interval.  Another  item  of  importance 
to  note  is  that  the  correlation  bandwidths  test  by  test  are  approximately 
equal  for  both  X-  and  C-bands .  This  observation  was  generally  true  for 
all  measurements  taken  at  all  the  overland  receiving  sites.  Point  Petre, 
the  over-water  site,  is  somewhat  different,  but  the  trend  (as  will  be 
shown  later)  is  present. 

Figures  65  and  66  are  representative  of  the  data  obtained  at  Whitford 
Field  and  were  taken  simultaneously  at  C-  and  X-bands.  Marked  breaks  in 
the  curves  did  not  occur  in  this  path  which  is  a  great  deal  longer  than 
the  Ontario  Center  path.  The  correlation  bandwidth  is  about  the  same  for 
either  frequency. 

Point  Petre,  the  overwater  path,  sometimes  deviates  completely  from 
what  mi^it  be  expected  in  that  Figure  67  shows  the  C-band  data  wider  than 
the  X-band  data  (Figure  68)  in  two  of  the  three  simultaneous  tests. 
Further,  a  study  of  the  Point  Petre  data  for  other  days  (Figures  69 
through  72)  show  essentially  the  trend  of  the  same  correlation  bandwidth 
for  each  frequency  band  with  the  slightly  wider  bandwidths  at  X-band. 

Correlation  bandwidth  appears  to  be  related  to  fade  rates  at  times, 
although  this  is  not  always  true.  For  example,  in  Figures  73  and  74,  it 
can  be  seen  that  the  correlation  bandwidth  test  that  has  the  highest  fade 
rate  (A)  nas  the  lowest  bandwidth  while  the  test  resulting  in  the  lowest 
fade  rate  (B)  has  the  widest  bandwidth.  This  trend  was  true  more  often 
than  not,  but  there  were  many  exceptions  as  illustrated  in  Figures 
75  and  76  in  which  approximately  the  reverse  is  observed.  Figure  77  is 
an  interesting  adjunct  to  Figure  63  for  studying  the  simultaneous  fade 
rates  versus  the  marked  discontinuity  in  correlation  bandwidth.  It  is 
noted  that  the  fade  rates  were,  in  general,  dropping  as  indicated  by 
tests.... 1,  2,  3,  and  4.  Between  tests  4  and  5  the  correlation  bandwidth 
reduced  to  about  half  its  previous  value  and  the  fade  rates  increased  for 
one  test  and  then  dropped  as  indicated  by  tests  ....5,  6,  and  8.  This 
observation  at  least  suggests  that  something  had  changed  within  the  common 
volume  and  adds  credence  to  the  proposal  that  the  atmosphere  in  the  common 
volume  had  changed  from  relatively  stable  to  unstable. 

The  tests  using  200  kHz  spacing  (Figures  78  and  79)  indicate  quite 
conclusively  that  the  derivative  of  the  correlation  coefficient  at  zero 
frequency  spacing  is  indeed  zero.  These  curves  are  typical  of  many  tests, 
all  of  which  yielded  essentially  the  same  results.  This  determination  is 
important  to  the  prediction  of  irreducible  error  rates  in  frequency  diver¬ 
sity  troposcatter  systems  which  according  to  Reference  21  results  in  many 
orders  of  magnitude  difference  in  predicted  performance. 
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Figure  65. 
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Figure  66. 
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Figure  68. 
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Figure  69. 
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Figure  70. 
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Figure  72. 
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Figure  74. 
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Figure  75. 
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Figure  76. 
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2. 


Overall  Correlation  Distributions 


Overall  correlation  coefficient  distributions  are  plotted  in  Figures 
80  through  95.  All  the  distribution  plots  for  correlation  coefficient 
data  are  included.  Only  wide  spaced  tests  are  used  in  these  plots.  Values 
of  the  envelope  cross  correlation  coefficient  are  shewn  on  the  left  axis. 
Corresponding  values  of  the  cross  correlation  coefficient,  p,  which  is  re¬ 
lated  to  the  envelope  cross  correlation  coefficient,  Pe,  by: 

Pe  =  P2 

21 

are  given  on  the  right  vertical  axis.  The  intersection  of  the  curves 
with  the  horizontal  line  at  Pe  =  0.4  defines  the  correlation  bandwidth. 

Both  the  upper  and  lower  deciles  as  well  as  the  median  correlation  for  the 
tests  included  are  shewn.  The  number  of  tests  in  each  distribution  is  also 
listed. 

Because  the  actual  antenna  beamwidth  w as  greater  for  C-band  (1.5  degrees 
HPBW)  than  X-band  (1.07  degrees  HPBW)  a  wider  correlation  bandwidth  would 
normally  be  expected  for  X-band.  However,  these  figures  illustrate  the 
general  finding  that  the  frequency  correlation  did  not  vary  markedly  from 
one  band  to  the  other.  In  fact,  for  the  Ontario  Center  summer  data,  the 
median  correlation  bandwidth  was  actually  slightly  larger  at  C-band. 

An  inverse  dependence  of  median  correlation  bandwidth  on  scatter 
angle  is  apparent  with  the  narrowest  medians  being  measured  for  Port  Byron 
and  Whitford  Field.  However,  no  systeratic  variation  of  the  medians  with 
time  of  year  is  evident  in  these  plots. 
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Figure  80. 
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Figure  81. 
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Figure  82. 
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Figure  83. 
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Figure  85. 
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Figure  87. 
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Figure  92. 
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Figure  93. 
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Figure  95. 
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3. 


Correlation  Bandwidth-Temperature  Effects 


In  a  format  similar  to  the  preceding  subsection.  Figures  96  through 
127  show  the  observed  correlation  function  distributions,  but  further 
subdivided  into  two  Fahrenheit  temperature  ranges.  These  temperature 
ranges  are  given  on  the  figures  and  were  chosen  so  that  the  crossover 
point  is  approximately  the  average  temperature  for  the  time  and  location 
of  the  measurements.  The  Whitford  Field  summer  X-band  data  included 
insufficient  tests  in  the  lower  temperature  range  (not  over  77  degrees) 
to  allow  a  determination  of  the  upper  and  lcwer  deciles. 

There  is  a  tendency  for  the  median  correlation  bandwidth  to  be  some¬ 
what  greater  for  the  lower  temperature  ranges  at  both  X-  and  C-band  and 
for  all  of  the  paths.  This  is  especially  noticeable  in  the  overwater  Point 
Petre  data  collected  during  the  September  testing  period. 
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Figure  98. 
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Figure  103. 
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Figure  108. 
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Figure  109. 
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Figure  113. 
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Figure  116. 
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Figure  117. 
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Figure  119. 
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Figure  124. 
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Figure  125. 
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4.  Correlation  Bandwidth-Diurnal  Effects 

In  a  format  similar  to  subsection  VI-A-2  ,  Figures  128  through  177 
give  the  observed  correlation  distributions,  but  subdivided  into  data  for 
four  6  hour  periods  through  the  day.  These  are  0001-0600,  0601-1200, 
1201-1800,  and  1801-2400  as  shown  in  the  figures.  For  several  of  the 
season-site- frequency  band-diumal  period  combinations  there  was  either  no 
or  insufficient  data  to  plot  medians.  For  others  there  was  too  little  data 
to  allow  a  plot  of  the  lower  and  upper  deciles. 

Hie  most  notable  feature  of  these  plots  is  the  lack  of  any  systematic 
variation  with  time  of  day.  Almost  any  stated  tendency  can  be  refuted  by 
several  counter  examples.  One  exception  lies  with  the  overwater  Point 
Petre  path  where  in  general,  at  both  X-  and  C-bands,  the  morning  or 
early  morning  correlation  bandwidths  were  slightly  greater  than  those 
observed  in  the  afternoon  or  evening. 
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Figure  139. 
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Figure  140. 
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Figure  145. 
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Figure  151. 
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Figure  152. 
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Figure  153. 
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Figure  156. 
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Figure  160. 
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Figure  161. 
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Figure  162. 
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Figure  163. 
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Figure  174. 
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Figure  176. 
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B.  FADE  RATE  DISTRIBUTIONS 


Figures  178  through  226  are  distributions  of  fade  rate  data.  Since 
more  rapid  fading  occurs  on  X-band  than  on  C-band ,  the  fade  rate  scale 
has  been  set  to  a  maximum  value  of  6C  Hertz  for  X-band  plots  and  20 
Hertz  for  C-band  plots.  The  ordinate  gives  the  percent  of  fades  less 
than  or  equal  to  the  fade  rate  given  by  the  abscissa.  Distribution  plots 
for  "All"  of  the  X-  and  C-band  tests  run  during  each  of  the  test  periods 
are  shown.  Typical  temperature  and  diurnal  effect  plots  are  also  included. 
When  no  marked  changes  occurred  with  time  and  temperature  variations,  the 
plots  au:e  not  included  in  this  report.  Furthermore,  after  several  examples 
of  typical  variations  of  tests  due  to  changing  conditions  are  given, 
additional  plots  are  omitted.  The  summary  bar  graphs  shown  at  the  end  of 
this  section  allow  a  comparison  of  the  relative  spread  of  all  plots. 

1.  Overall  Fade  Rate  Distributions 

Figures  178  through  193  give  the  distributions  of  fade  rates  observed 
over  the  four  paths.  The  plots  are  arranged  in  the  order  of  season  or 
month,  then  by  site,  and  finally  by  frequency  band.  For  all  curves  the 
observed  upper  and  lower  deciles  and  median  are  indicated  together  with 
the  number  of  tests  for  which  the  distributions  were  calculated. 

For  all  paths  the  fade  rate  distributions  are  very  similar  and  indi¬ 
cate  no  seasonal  dependence.  The  variation  in  the  deciles  is  not  statis¬ 
tically  significant.  There  is  however  a  marked  difference  between  the  X- 
and  C-band  data.  On  the  average  50  percent  of  all  X-band  fades  occurred 
at  a  rate  less  than  about  4  Hz  while  at  C-band  the  value  is  approximately 
3  Hz.  At  the  80  percent  level  the  values  are  approximately  9  Hz  for  X- 
band  and  6  Hz  for  C-band. 
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Figure  184. 
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Fade  Rate  Distributions  -  Temperature  Effects 


In  a  format  similar  to  the  preceding  subsection,  Figures  194  through 
209  show  the  observed  fade  rate  distributions,  but  subdivided  into  two 
Fahrenheit  temperature  ranges.  These  temperature  ranges  jure  given  on  the 
figures  and  were  chosen  so  that  the  crossover  point  is  approximately  the 
average  temperature  for  the  time  and  location  of  the  measurements. 

The  first  six  figures  were  calculated  from  X-band  data,  the  last  ten 
from  C-band.  Only  a  representative  sample  of  available  plots  is  presented 
here  since  all  data  collected  indicated  there  was  no  significant  dependence 
of  fade  rate  on  surface  temperature. 
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Figure  195. 


150 


8 
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Figure  202. 
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3. 


Fade  Rate  Distributions  -  Diurnal  Effects 


In  a  format  similar  to  subsection  VI-B-1,  Figures  210  through  226 
give  the  observed  fade  rate  distributions  but  subdivided  into  data  for 
four  six  hour  periods  through  the  day.  These  are  0001-0600,  0601-1200, 
1201-1800,  and  1801-2400  as  shown  in  the  figures.  The  data  included  in 
three  of  the  figures  contained  too  few  tests  to  allow  a  computation  of 
the  deciles. 

The  first  three  figures  are  based  on  X-band  data,  the  last  14  on 
C-band.  Only  a  representative  sample  is  given  here.  A  common  feature 
of  many  of  the  plots,  including  those  not  shown  here,  is  a  tendency  for 
the  median  fade  rate  at  both  X-  and  C-band  to  be  slightly  greater  during 
the  evening  hours.  There  are  exceptions  to  this  finding,  however,  as  is 
evident  in  the  Point  Petre  data. 
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Figure  213. 
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Figure  216. 
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Figure  219. 
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Figure  220. 
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Figure  221. 
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C.  DEPTH  OF  FADE  DISTRIBUTIONS 

Fade  depth  distributions  are  shewn  in  Figures  227  through  253.  On 
all  of  these  plots  cumulative  percentage  is  plotted  on  a  normal  proba¬ 
bility  scale.  A  fade  was  indicated  by  any  drop  in  level  from  the  median 
signal  level,  the  median  being  established  for  blocks  of  2000  samples. 
Only  typical  examples  of  temperature  and  diurnal  effects  on  these  kinds 
of  data  are  shown.  For  those  plots  which  are  omitted  from  this  report, 
either  the  spread  of  values  corresponds  closely  to  the  included  "All" 
plots  or  the  range  of  values  is  similar  to  the  included  temperature 
and  diurnal  effect  plots.  The  comparative  spread  of  values  is  shown  on 
the  summary  bar  graphs  at  the  end  of  this  section. 

1.  Overall  Depth  of  Fade  Distributions 

Figures  227  through  242  give  the  overall  fade  depth  distributions 
over  the  four  paths.  In  addition  to  the  median  for  all  tests  included, 
the  upper  and  lower  deciles  are  plotted  where  justified  by  sufficient 
data.  The  number  of  tests  for  each  distribution  is  shown  on  the  figures. 

For  all  paths  the  depth  of  fade  distributions  are  similar  and  show 
no  seasonal  or  frequency  dependence.  The  observed  median  fade  depth 
was  approximately  11  dB.  On  the  average  less  than  10  percent  of  all 
fades  exceeded  20  dB. 
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Figure  235. 
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Figure  241. 
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Figure  242. 
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2.  Depth  of  Fade  Distributions -Temperature  Effects 

In  a  format  similar  to  the  preceding  subsection ,  Figures  243  through 
246  show  the  observed  depth  of  fade  distributions  but  subdivided  into 
two  Fahrenheit  temperature  ranges.  These  temperature  ranges  are  given  on 
the  figures  and  were  chosen  so  that  the  crossover  point  is  approximately 
the  average  temperature  for  the  time  and  location  of  the  measurements. 

The  figures  shown  are  representative  and  indicate  no  temperature 
effect  on  the  fade  depths. 
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Figure  243. 
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Depth  of  Fade  Distributions  -  Diurnal  Effects 


In  a  format  similar  to  that  of  subsection  VI-C-1,  Figures  247  through 
253  give  the  observed  depth  of  fade  distributions  but  subdivided  into  data 
for  four  6  hour  periods  through  the  day.  These  are  0001-0600,  0601-1200, 
1201-1800,  and  1801-2400  as  shewn  in  the  figures.  The  data  included  in 
four  of  the  figures  contained  too  few  tests  to  allow  a  computation  of  the 
deciles. 

Only  a  representative  sample  of  all  of  the  available  data  is  pre¬ 
sented  here  since  there  was  no  indication  of  any  diurnal  effect. 
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Figure  251. 
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Figure  253 


D.  FADE  DURATION  DISTRIBUTIONS 


Fade  duration  distributions  are  shown  in  Figures  254  through  286  for 
fade  levels  of  0,  5,  10,  15,  and  20  dB  below  the  median.  On  all  of  these 
plots  cumulative  percentage  is  plotted  on  a  normal  probability  scale.  Only 
typical  examples  of  temperature  and  diurnal  effects  on  these  kinds  of  data 
are  shown.  For  those  plots  which  are  omitted  from  this  report,  either  the 
spread  of  values  corresponds  closely  to  the  "All"  plots  which  are  included 
or  the  range  of  values  is  similar  to  the  temperature  and  diurnal  effect 
plots  which  are  included.  The  comparative  spread  of  values  is  shown  on  the 
summary  bar  graphs  at  the  end  of  this  section. 

1.  Overall  Fade  Duration  Distributions 

Figures  254  through  269  give  the  overall  fade  duration  distributions 
over  the  four  paths.  The  plots  are  arranged  in  the  order  of  season  or 
month,  then  by  site,  and  finally  by  frequency  band.  The  number  of  tests 
for  each  distribution  is  shown  on  the  figures.  The  fade  duration  dis¬ 
tribution  is  given  for  the  median  and  four  other  levels. 

For  all  paths  the  fade  duration  distributions  are  similar  and  no 
seasonal  dependence  is  evident.  As  might  be  expected  from  the  fade  rate 
distributions,  there  is  a  difference  between  the  X-  and  C-band  durations. 
For  the  fade  duration  at  the  median  level  the  average  value  at  X-band  is 
approximately  80  milliseconds  while  at  C-band  the  corresponding  figure  is 
110  milliseconds. 
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Figure  254. 


T 


DISTRIBUTION  OF  FADE  DURATION 

ONTARIO  CENTER,  SUMMER,  C-BflNO 

RLU-46  TESTS 

Figure  255. 
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Figure  256. 
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Figure  257 
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Figure  259. 
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Figure  261. 
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Figure  264. 
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Figure  265 
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Figure  266. 
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Figure  267. 
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Figure  268. 
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Figure  269. 
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2. 


Fade  Duration  Distributions  -  Temperature  Effects 


In  a  format  similar  to  the  preceding  subsection,  Figures  270  through 
275  show  the  observed  depth  of  fade  distributions  but  subdivided  into  two 
Fahrenheit  temperature  ranges.  These  tomperature  ranges  are  given  on  the 
figures  and  were  chosen  so  that  the  crossover  point  is  approximately  the 
average  temperature  for  the  time  and  location  of  the  measurements. 

The  figures  shown  are  representative  and  indicate  no  temperature 
effect  on  fade  duration. 
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Figure  270. 
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Figure  271. 
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3.  Fade  Duration  Distributions  -  Diurnal  Effects 

In  a  format  similar  to  that  of  subsection  VI-D-1,  Figures  276 
through  286  give  the  observed  fade  duration  distributions  but  subdivided 
into  data  for  four  6  hour  periods  through  the  day.  These  are  0001-0600, 
0601-1200,  1201-1800,  and  1801-2400  as  shown  in  the  figures. 

Only  a  representative  sample  of  all  of  the  available  data  is  pre¬ 
sented  here  sinc<3  there  was  no  indication  of  any  diurnal  effect. 
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Figure  277. 
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Figure  278. 
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Figure  279. 
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Figure  281. 
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Figure  282. 


DISTRIBUTION  OF  FAOE  DURATION 

ONTARIO  CENTER.  XTOBER.  C-BANO 
FROM  1801  TO  2400.  25  TESTS 


Figure  283. 
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E.  SIGNAL  AMPLITUDE  DISTRIBUTIONS 

Signal  amplitude  distributions  are  shown  in  Figures  287  through  324. 
On  all  of  these  plots,  cumulative  percentage  is  plotted  on  a  normal  proba¬ 
bility  scale  against  variations  in  signal  strength  from  the  median.  Only 
typical  examples  of  temperature  and  diurnal  effects  on  these  kinds  of 
data  are  shewn.  For  those  plots  which  are  omitted  from  this  report, 
either  the  spread  of  values  corresponds  closely  to  the  "All"  plots  which 
are  included  or  the  range  of  values  is  similar  to  the  temperature  and 
diurnal  effect  plots  which  are  included.  The  comparative  spread  of 
values  is  shown  on  the  summary  bar  graphs  at  the  end  of  this  section. 

1.  Overall  Signal  Amplitude  Distributions 

Figures  287  through  302  give  the  overall  signal  amplitude  distri¬ 
butions  over  the  four  paths.  The  plots  are  arranged  in  the  order  of 
season  or  month,  then  by  site,  and  finally  by  frequency  band.  The  number 
of  tests  for  each  distribution  is  shown  on  the  figures.  For  all  cases 
where  sufficient  data  was  available  the  upper  and  lower  deciles  have  been 
plotted  in  addition  to  the  medians. 

For  all  paths  the  signal  amplitude  distributions  are  similar  and  no 
seasonal  dependence  is  evident.  The  median  curves  as  plotted  on  this 
log-normal  scale  yield  distributions  that  are  usually  Rayleigh. 
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Figure  288. 
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Figure  290. 
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Figure  294. 
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Figure  296. 
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Figure  300. 
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Figure  302. 
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2.  Signal  Amplitude  Distributions  -  Temperature  Effects 


In  a  format  similar  to  the  preceding  subsection.  Figures  303  through 
310  show  the  observed  signal  amplitude  distributions  but  subdivided  into 
two  Fahrenheit  temperature  ranges.  These  temperature  ranges  are  given 
on  the  figures  and  were  chosen  so  that  the  crossover  point  is  approximately 
the  average  temperature  for  the  time  and  location  of  the  measurements. 

The  figures  shown  are  representative  and  indicate  no  temperature  effect 
on  the  signal  amplitude  distribution. 
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Figure  304. 
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Figure  306. 
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Figure  310. 
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3.  Signal  Amplitude  Distributions  -  Diurnal  Effects 

In  a  format  similar  to  that  of  subsection  Vl-E-1,  Figures  311  through 
324  shew  the  observed  signal  amplitude  distributions  but  subdivided  into 
four  6  hour  periods  through  the  day.  These  are  0001-0600,  0601-1200, 
1201-1800,  and  1801-2400  as  shown  in  the  figures. 

Only  a  representative  sample  of  all  of  the  available  data  is  presented 
here  since  there  was  no  indication  of  any  diurnal  effect.  The  ujiper  decile 
Ontario  Center,  October,  C-band  (Figure  318)  is  not  statistically  significant 
since  it  was  contributed  by  only  two  tests. 
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Figure  314. 
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Figure  318. 
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Figure  322. 
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F.  SUMMARY  BAR  GRAPHS 


Figures  325  through  360  are  summary  bar  graphs  for  X-  and  C-band 
data  from  the  four  sites.  Points  read  from  the  distribution  plots  are 
presented  in  this  form  to  permit  comparison  of  the  range  of  values  during 
different  test  conditions.  Where  sufficient  data  exists,  the  upper  and 
lower  deciles  are  plotted  in  addition  to  the  median.  Summary  bar  graphs 
are  not  shown  for  fade  depth  or  duration  or  signal  amplitude,  for  Port 
Byron  X-band  due  to  lack  of  sufficient  fade  margin. 

With  the  exception  of  the  fade  duration  summaries,  the  figures  give 
values  for  the  upper  and  lower  deciles  as  well  as  the  median.  In  cases  of 
insufficient  data  only  the  medians  are  shown.  Irrespective  of  the  type  of 
chart,  the  median  is  to  be  interpreted  as  the  median  value  for  all  included 
tests  for  the  particular  parameter  displayed.  The  deciles  are  contributed 
by  the  averages  in  the  upper  and  lower  10  percent  of  all  included  tests 
which  individually  were  in  the  upper  or  lower  10  percent  of  the  parameter 
displayed. 

Figures  325  through  331  give  the  summaries  of  the  observed  frequency 
correlation  coefficients.  The  ordinate  gives  the  correlation  bandwidth, 
defined  by  =  0.4,  in  MHz.  As  an  example,  in  Figure  325  the  median 
correlation  bandwidth  at  Ontario  Center  in  the  summer  for  all  X-band  data 
is  seen  to  be  approximately  3  MHz.  That  is,  each  of  the  81  tests  included 
is  characterized  by  a  value  of  Pe  and  the  median  of  these  is  3  MHz.  The 
average  of  the  8  tests  (10%)  having  the  highest  individual  values  of  Pe 
is  the  upper  decile  value,  here  about  8.8  MHz.  Likewise,  the  lower  decile 
is  about  1.5  MHz. 

Figures  332  through  339  give  the  fade  rate  summaries.  The  ordinate 
gives  the  observed  median  fade  rate  in  Hz.,  i.e.,  the  50%  level.  The 
plotted  medians  are  actually  medians  of  medians.  For  example,  in  Figure 
332  for  Ontario  Center  in  the  summer  for  all  X-band  data  50%  of  all  fades 
were  less  than  2.4  Hz.  Thus  of  the  88  tests  comprising  this  summary  each 
had  a  median  fade  rate.  The  median  of  these  is  2.4  Hz.  The  upper  and  lower 
deciles  are  given  by  the  averages  in  the  upper  and  lower  10  percent  of  tests 
(each  9  tests)  each  individually  had  the  highest  and  lowest  median  (50%) 
fade  rates. 

Figures  340  through  346  give  the  fade  depth  summaries.  The  ordinate  is 
the  median  fade  depth  in  dB,  i.e.,  the  50%  level.  As  for  the  fade  rate 
summaries  described  above,  the  plotted  medians  are  actually  medians  of 
medians  and  a  similar  interpretation  to  that  given  above  applies  to  the 
deciles. 

Figures  347  through  353  give  the  fade  duration  summaries.  The  ordinate 
is  the  median  fade  duration  in  milliseconds  for  all  included  tests  and 
the  three  levels  shown  (0,  10,  and  20  dB)  are  with  respect  to  the  individual 
median  signal  amplitude  levels.  As  an  example,  again  choosing  the  Ontario 
Center,  summer,  X-band,  all  data  shown  in  Figure  347,  the  median  fade  dura¬ 
tion  below  the  individually  measured  median  signal  amplitudes  is  approxi¬ 
mately  105  milliseconds.  At  levels  10  dB  below  the  individually  measured 
median  signal  amplitudes,  the  median  fade  duration  is  35  milliseconds. 
Similarly,  at  the  20  dB  level  the  median  fade  duration  is  about  15  milli¬ 
seconds  . 

Finally,  the  signal  amplitude  summaries  are  shown  in  Figures  354 
through  360.  The  ordinate  gives  the  signal  amplitude  value,  in  dB  below 
the  individual  test  median  signal  amplitudes,  which  was  exceeded  95  percent 
of  the  test  time.  The  plotted  upper  decile  indicates  that  the  upper  10 
percent  of  the  tests  included  had  signal  amplitudes  which  varied  from  the 
individual  median  amplitudes  by  no  more  than  the  ordinate  value  for  95  percent 
of  the  time.  Similarly,  the  10  percent  of  included  tests  which  displayed 
relatively  deeper  fading  comprise  the  lower  decile.  The  individual  median 
signal  amplitudes  on  an  absolute  scale  naturally  varied  from  test  to  test. 
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Figure  327. 
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Figure  328. 
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Figure  329. 
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Figure  330. 
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FADE  DEPTH  SUMMARY 
Figure  342. 
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FflOE  DEPTH  SUMMRRY 
Figure  344. 
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FROE  DEPTH  SUM 
Figure  346 


249 


FflOe  DURATION 
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Figure  350. 
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Figure  352. 


FfiCE  DURATION  SUMMARY 
Figure  353. 


SIGNAL  AMPLITUDE  SUMMARY 
Figure  356. 
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SIGNAL  AMPLITUDE  SUMMARY 
Figure  358. 
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SIGNAL  AMPLITUDE  SUMMARY 


SIGNAL  AMPLITUDE  SUMMARY 
Figure  360. 


G.  SIGNIFICANT  FIELD  OBSERVATIONS  AND  UNUSUAL  PHENOMENA 

Some  unusual  phenomena  observed  in  the  field  tests  were  ducting  and 
its  associated  high  signal  strengths.  Ducting  on  C-band  and  a  simultaneous 
normal  situation  on  X-band  was  noticed  on  one  occasion.  High  fade  rates 
were  observed  many  times  due  to  propagation  phenomena  and  due  to  aircraft. 
Occasionally  there  were  unusual  shapes  in  the  correlation  coefficient 
curves  that  suggest  a  (sin  x)/x  function  versus  frequency.  There  appear 
to  be  relationships  between  the  cloud  ceiling  height  and  the  received 
signal  strengths  and  also  the  cloud  ceiling  height  and  the  correlation 
bandwidths  measured. 

1 .  Ducting 

Ducting  signals  are  identified  by  their  relatively  high  signal 
strengths  and  often  low  depth  of  fades  accompanied  usually  by  a  wide 
correlation  bandwidth.  However,  the  high  signal  strength  is  the  only 
common  denominator  observed  in  all  ducting  situations  identified  during 
these  tests.  Figures  361  through  364  show  a  case  of  X-band  ducting 
observed  at  Ontario  Center  on  11  August  at  2035.  (The  C-band  transmitter 
was  not  operating  at  this  time.)  The  correlation  bandwidth  is  seen  to 
be  very  wide,  but  the  fades  on  which  the  computation  was  made  are  very 
shallow  and  seldom  more  than  5  dB  below  the  median.  The  fade  rate  is 
enormous  with  10  percent  of  the  fades  greater  than  60  Hz.  In  this  situation 
the  signal  strengths  were  in  the  vicinity  of  -60  dBm.  At  2300  on  this 
same  date  ducting  occurred  again  with  the  C-band  results  shown  in  Figures 
365  through  368.  Here  the  correlation  bandwidth  is  narrow  with  some  fades 
greater  than  5  dB  below  the  median.  The  signal  strength  was  about  -66  dBm. 

Usually  ducting  was  observed  on  the  Ontario  Center  and  the  Point  Petre 
paths  in  the  stammer  and  early  fall.  This  phenomena  was  never  identified 
on  the  Whitford  Field  or  Port  Byron  paths.  At  Whitford,  the  antenna  take¬ 
off  angle  was  probably  too  large  for  ducting  to  occur,  while  on  the  Port 
Byron  path  the  winter  conditions  and  short  testing  period  reduced  the 
probability  of  observing  any  ducting  phenomena.  Ducting  usually  occurred 
at  night  from  about  two  hours  after  sunset  until  about  an  hour  before 
sunrise.  In  all  cases  the  C-band  went  into  ducting  prior  to  X-band  and 
remained  after  the  X-band  signal  had  returned  to  normal.  In  one  instance 
over  the  Point  Petre  path,  a  ducting  condition  produced  a  signal  strength 
of  -50  dBm  on  C-band  with  a  transmitted  power  of  only  6  watts. 

Another  case  of  ducting  occurred  at  Point  Petre  where  evidently  the 
C-band  was  experiencing  a  more  classic  case  while  X-band  was  not  fully 
ducting.  Figures  369  through  373  show  the  C-band  propagation  and  Figures 
374  through  378  show  the  X-band.  The  signal  strength  was  -76  dBm  for 
X-band  and  -66  dBm  for  C-band.  The  C-band  ducting  had  the  wide  correlation 
bandwidth,  very  high  fade  rates  with  mostly  shallow  fades  of  which  80  per¬ 
cent  were  less  than  4  dB.  The  unexpected  phenomenon  is  that  the  X-band 
was  not  strictly  ducting.  True,  it  had  a  wide  correlation  bandwidth  most 
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of  the  time/  but  the  fade  rates  were  within  typical  bounds  and  so  were 
the  signal  amplitude,  fade  duration,  and  depth  distributions.  In  this 
type  of  environment  the  signal  is  of  such  strength  that  both  frequency¬ 
time  and  adoptive  frequency  modems  operate  very  well.  The  loss  of  diver¬ 
sity  gain  is  not  important. 
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CORRELATION  COEFFICIENT  (RHO  El/2) 


Figure  364. 
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ENVELOPE  CROSS  CORRELATION  COEFFICIENTS 

ONTARIO  CENTER*  SUMER 

C-BRND-N10E  DUCTING  Figure  365. 


FADE  RATE  DISTRIBUTION 

ONTARIO  CENTER.  UtO 

C-flAND  DUCTING  Figure  366. 
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PERCENT  (P  FROES  fts?THAlFoR  E$JflP  TO  FRo£  DURATION 
DISTRIBUTION  OF  FADE  DURATION 

ONTARIO  CENTER*  SUMMER  0611290011 

C-6RND  DUCTING  Figure  367. 
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SIGNAL  AMPLITUDE  LEVEL 

ONTRRIO  CENTER.  SUMMER 

C-flRND  OUCTING  Figure  368. 
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ENVELOPE  CROSS  CORRELATION  COEFFICIENTS 


7-SO  10-00  12-60 

FADE  RATE  IN  HERTZ 
FADE  RATE  DISTRIBUTION 


POINT  PETRE*  SEPT 

C-*ND  Figure  370. 
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SIGNAL  STRENGTH ,  DB  ABOVE  MEDIAN  .  Ffloe  DURATION  IN  MILLISECOfOS 
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DISTRIBUTION  OF  FADE  DURATION 

POINT  PETRE,  8EPT 

C-flfH)  0923-4  TO  8  Figure  371. 
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POINT  PETRE,  8EPT 

C-flRND  Figure  372. 
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ENVELOPE  CROSS  CORRELATION  COEFFICIENT  (RHO  E)  DEPTH  OF  FADE »  DB 
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CORRELATION  COEFFICIENT 


1 


DISTRIBUTION  OF  FADE  DURATION 


POINT  PETRE.  8EPT 


X-flflND  0983-4  TO  8  Figure  376. 
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X-BflND  Figure  378. 
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2.  Aircraft  Effects 

Aircraft  effects  should  not  be  called  an  anomalous  condition  _ 

it  is  a  condition  often  present  in  both  military  and  civilian  environments. 
The  presence  of  the  aircraft  for  a  short  interval  has  very  little  effect 
on  a  statistical  test  of  relatively  long  time  duration.  Figures  379  and 
380  show  the  presence  of  the  aircraft  by  the  rise  in  correlation  coeffi¬ 
cient,  but  has  little  effect  on  the  overall  test.  The  fade  rates  shown 
in  Figures  381  and  382  are  averaged  over  the  entire  test  number  15  and  13 
respectively  to  indicate  that  the  effect  is  noticeable  over  the  entire 
test.  During  the  time  that  the  aircraft  were  actually  present  in  the  data, 
fade  rates  were  much  higher  than  shown  on  the  curves. 

The  frequency- time  modems  can  operate  through  this  type  of  environment. 
The  adaptive  frequency  modems  were  previously  found  to  have  difficulty  due 
to  the  high  fade  rates  (Reference  21)  . 


FREQUENCY  SEPARATION  IN  MHZ 


ENVELOPE  CR08S  CORRELATION  COEFFICIENTS 

GNTHtI0  CCMTBt»  ftJMCX 
XHRO'NIK  AINUNE  EFFECT  I 

Figure  379. 


COEFFICIENT  ( 


The  signatures  of  the  aircraft  presence  are  of  interest  because  the 
resulting  simultaneous  C-  and  X-band  AGC  versus  time  curves  are  decidedly 
different.  Figure  383  illustrates  this  phenomena  clearly.  The  C-band 
signals  receive  interference  first  and  this  interference  remains  after 
the  X-band  interference  has  disappeared.  This  observation  appears  to  be 
due  to  the  combined  effects  of  the  wider  antenna  beamwidth  and  longer 
wavelength  of  the  C-band  signals  which  causes  the  aircraft  to  become 
illuminated  prior  to  the  X-band  illumination  with  fades  proportional  to 
the  wavelengths  involved  relative  to  the  aircraft  motion  involved. 


ENVELOPE  CROSS  CORRELATION  COEFFICIENTS 
ONTPRIO  CENTER.  HNO 
X-aM'NlDE  RUftftC  EFFECT  2 


Figure  380. 
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C-band 


Figure  383.  Simultaneous  AGC  Recording  of  Aircraft  Effects  in 

a  Tropo  Path  During  Test  02101430W 


3.  Unusual  Shape  in  the  Correlation  Coefficient  Curves 

Some  unusual  shapes  in  the  correlation  coefficient  curves  were  noted 
often  enough  to  consider  them  as  perhaps  a  part  of  the  fading  and  multi- 
path  mechanism  occasionally  present  in  the  common  volume.  This  phenomenon 
was  not  observed  in  the  previous  troposcatter  propagation  work  in  Ref¬ 
erence  21.  If  it  were  present,  it  went  undetected.  The  correlation  band¬ 
width  curves  are  presented  in  Figures  384  and  385  because  they  might 
give  a  clue  to  the  behavior  of  the  common  volume.  The  curves  appear  as  if 
they  are  following  some  sort  of  (sin  x)/x  function.  In  Figure  384,  it 
appears  suddenly  in  test  12  at  1930  hours  and  gradually  diminishes, 
changing  its  periodicity,  over  the  several  hours.  The  phenomenon  is 
repeated  to  a  lesser  degree  in  Figure  385. 
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ENVELOPE  CROSS  CORRELATION  COEFFICIENT  (RHO  E)  ENVELOPE  CROSS  CORRELATION  COEFFICIENT  (RHO  E) 


CORRELATION  COEFFICIENT  (RHO  El/2)  CROSS  CORRELATION  COEFFICIENT  (RHO  El/2) 


PERCENT 


4.  High  Fade  Rates 

Occasionally  very  high  fade  rates  were  observed.  As  discussed  in  a 
subsequent  section,  these  could  have  some  effect  on  adaptive  frequency 
modems.  Figures  386  and  387  are  examples  of  some  of  the  very  high  fade 
rates.  These  high  fade  rates  render  the  best  frequency  select  circuits  to 
be  spoofed  and  to  cause  the  adaptive  frequency  modem  to  go  to  a  less  than 
optimum  frequency.  Sometimes  this  can  cause  an  error  in  the  received  fre¬ 
quency  commands  and  can  cause  the  link  to  become  broken. 


5-00  7*50  10-00  12-50 

FADE  RATE  IN  HERTZ 

FADE  RATE  DISTRIBUTION 
0NTRRI0  CENTER »  8UWER 
C-BflND 

Figure  386. 


279 


u-00 


7.80  lkoo j£-60 30-00 37-50 

FADE  RATE  IN  HERTZ 
FADE  RATE  DISTRIBUTION 


ONTARIO  CENTEX.  BUHNER 
X-flRND 


46-00 


5^ -BO  60-00 


5.  Effects  of  Cloud  Ceiling  on  Troposcatter 

It  was  observed  in  the  field  that  the  ceiling  seemed  to  affect  the 
received  signal  strength.  Weather  data  was  examined  to  find  days  with 
all  meteorological  data  remaining  relatively  constant  except  the  ceiling. 
TWo  such  days  were  found  and  plots  of  signal  strength  versus  ceiling  and 
correlation  bandwidth  versus  ceiling  are  shewn  in  Figures  388  and  389, 
respectively.  The  signal  strength  versus  ceiling  data  indicated  that 
as  the  ceiling  height  increased  the  signal  strength  increased  until  the 
ceiling  was  equal  to  or  slightly  greater  than  the  intersection  of  the  lower 
rays  and  then  began  to  decrease  as  the  ceiling  height  continued  to  rise. 

The  correlation  bandwidth  versus  ceiling  data  indicated  that  the  corre¬ 
lation  bandwidth  was  wide  when  the  ceiling  was  low  and  narrow  when  the 
ceiling  was  high.  The  plot  was  similar  to  a  hyperbolic  function  with 
the  knee  of  the  curve  occurring  at  the  altitude  of  the  intersection  of  the 

lower  rays. 

These  observations  are  by  no  means  to  be  considered  as  conclusive  at 
this  time,  for  the  quantity  of  data  used  in  the  observations  is  insuf¬ 
ficient.  They  are,  however,  very  interesting  observations  and  might  be 
considered  for  further  study  in  the  future. 


> 


280 


Correlation  1/2  Bandwidth,  MHz 


Ceiling  Altitude,  Thousands  Of  Feet 
Figure  388.  Correlation  Bandwidth  versus  Ceiling  Height 


Figure  389.  Signal  Strength  versus  Ceiling 
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Envelope  Correlation  1/2  Bandwidth  MHz  p  =  0.4 


Figures  390  through  406  summarize  the  troposcatter  characteristics  observed 
over  the  Tobyhanna,  Pennsylanvia  to  Hexagon,  Fort  Monmouth,  New  Jersey  link 
in  the  period  February  through  July,  1968.  These  tests  were  made  at  C-band 
and  employed,  for  most  of  the  test  period,  10  foot  parabolic  antennas  at 
both  the  transmitter  and  receiver.  Full  details  of  the  link  and  the  obser¬ 
vations  are  given  by  Branham  et  al.  in  "Troposcatter  Transmission  of  High 
Speed  Digital  Signals"  (reference  21).  In  brief,  this  link  may  be  described 
as  a  94  mile  path  over  predominantly  mountainous  terrain  which  becomes 
coastal  near  the  New  Jersey  termination.  The  scatter  angle  is  0.674  degrees 
which  is  smaller  than  that  of  any  other  path  tested  in  this  program. 

The  figures  follow  the  same  nomenclature  used  in  presenting  the  pre¬ 
ceding  results  and  are  included  primarily  for  reference  purposes.  It  is 
seen  that  the  median  correlation  bandwidth  for  all  data  is  approximately 
4.4  MHz  (Figure  390).  This  is  approximately  1  to  2  MHz  wider  than  the  med¬ 
ian  correlation  bandwidths  observed  in  the  N.  Y.  paths.  The  wide  correla¬ 
tion  bandwidth  implies  a  smaller  multipath  spread  and  therefore  hints  that 
the  effective  beamwidth  over  this  link  was  still  less  than  the  actual  (the 
model  of  Section  II  using  equation  43  would  set  the  effective  beamwidth 
equal  to  the  actual  beamwidth  for  this  scatter  angle) .  The  fade  rates 
(Figure  404)  are,  however,  comparable  to  those  observed  here. 


H.  DATA  FROM  TOBYHANNA/HEXAGON  PATH  (PREVIOUS  PROGRAM) 
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Figure  390.  Seasonal  and  Diurnal  Effects  on  Cross-correlation  Coefficient 
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Frequency  Separation  (A  f)  in  MHz 

Figure  391.  Envelope  Cross-correlation  Coefficients,  Entire  Program  (ALL)  682  Tests 
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Figure  394.  Envelope  Cross -correlation  Coefficients,  Winter  (ALL)  79  Tests 
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Frequency  Separation  (A  f)  in  MHz 


2/ja 0  juaxotjjaoo  uoijBiajaoo  ssoao 


o  ®  oo  r-  co  m  Tf  oo  n  © 


ao 


r- 


(O 


m 


00 


M 


O 


288 


Frequency  Separation  (A  f)  in  MHz 

Figure  396.  Envelope  Cross-correlation  Coefficients,  Winter  (PM)  34  Tests 


6'0 


289 


Figure  397.  Envelope  Cross -correlation  Coefficients,  Spring  (ALL)  328  Tests 


Frequency  Separation  (A  f)  in  MHz 

Figure  398.  Envelope  Cross -correlation  Coefficients,  Spring  (AM)  142  Tests 
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Frequency  Separation  (Af)  in  MHz 

Figure  399.  Envelope  Cross -correlation  Coefficients,  Spring  (PM)  186  Tests 
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Frequency  Separation  (A  f)  in  MHz 

Figure  402.  Envelope  Cross -correlation  Coefficients,  Summer  (PM)  165  Tests 
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Figure  403.  Fade  Duration  Distribution,  Entire,  ALL 
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Figure  404.  Fade  Rate  Summary 
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Figure  405.  Fade  Rate  Distribution,  Summer,  AM 
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Figure  406.  Signal  Amplitude  Distribution,  Entire,  ALL> 


I.  EXPECTED  PERFORMANCE  OF  VARIOUS  MODEMS 

In  a  previous  MALLARD  program  several  frequency  diversity  modems 
were  evaluated  over  the  Tobyhanna/Hexagon  troposcatte-:  link.  These  were 
reported  in  Reference  21.  The  diversity  obtained  by  the  use  of  frequency 
separation  between  branches  depended  on  decorrelated  fading  between  them. 

On  the  other  hand,  the  data  handling  capacity  of  any  individual  branch 
depended  on  the  correlation  bandwidth  of  the  branch.  Thus  it  was  found  that 
the  very  mechanism  that  provided  frequency  diversity  also  was  the  mechanism 
that  limited  the  channel  capacity  due  to  multipath  distortions  brought  on 
by  propagation  through  dispersive  media.  Nevertheless,  it  was  found  that 
frequency  diversity  at  data  rates  of  slightly  over  one  megabit  was  obtain¬ 
able  through  the  use  of  simple  frequency-time  codes.  An  adaptive  frequency 
selection  system  was  also  tested.  It  employed  one  branch  at  a  time  with  a 
selection  of  frequency  on  a  best-by-test  basis  and  a  transfer  to  the  apparent 
best  frequency  through  the  use  of  frequency  change  command  codes  transmitted 
across  the  link. 

From  the  data  obtained  during  this  test  program  it  appears  that  the 
order  of  frequency  diversity  at  X-band  is  essentially  the  same  as  that 
obtainable  at  C-band  over  these  paths  with  10  foot  antennas.  If  the 
antenna  sizes  were  to  be  increased  one  might  expect  the  correction  band¬ 
width  of  X-band  to  eventually  increase  as  the  beamwidth  became  less  than 
maximum  usable  effective  beamwidth.  Eventually,  further  increases  in 
antenna  size  would  cause  the  C-band  correlation  bandwidth  to  also  increase 
for  the  same  reason.  Increases  in  correlation  bandwidth  would  increase 
the  data  rates  per  branch  of  frequency  diversity,  but  reduce  the  diver¬ 
sity  obtainable  among  the  individual  branches. 

The  frequency- time  modems  can  be  expected  to  perform  as  well  at 
X-band  as  they  do  at  C-band  from  the  standpoint  of  diversity  obtained. 

Some  fade  margin,  however,  is  lost  at  X-band  due  to  increased  path  loss 
and  medium  to  aperture  coupling  losses  at  the  higher  frequency.  The  high 
fade  r?tes  obtained  over  these  paths  at  both  frequencies  do  not  appear  to 
pose  a  problem  for  the  frequency- time  modems  as  they  can  have  high  speed 
AGC  circuits  that  can  cope  with  the  rapid  changes  in  signal  strength. 
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The  adaptive  frequency  modems  should  obtain  the  same  orders  of 
diversity  and  path  losses  as  the  frequency-time  modems.  They  do,  however, 
have  a  problem  in  coping  with  enormous  fade  rates  that  have  been  encoun¬ 
tered  over  these  paths.  There  is  a  definite  time  required  between  the 
tests  to  determine  the  apparent  best  frequency  and  the  effecting  of  a 
change  to  that  frequency.  This  requires  the  propagation  across  the  link 
of  the  best  frequency  command  code  twice  (minimum)  so  that  even  before 
the  system  can  transfer  to  the  new  best  frequency  it  may  no  longer  be 
optimum.  The  system  tested  during  the  previous  program  must  therefore 
be  further  developed  to  operate  through  higher  fade  rates. 

All  modems  should  operate  through  a  ducting  condition  since  it  is 
characterized  by  an  unusually  high  signal  strength  with  shallow  fades. 
Frequency  diversity  is  not  readily  obtainable  during  ducting,  but  with 
such  strong  signals  the  diversity  gain  is  usually  not  needed  anyway,  and 
all  modems  operate  essentially  error  free.  The  apparent  small  corre¬ 
lation  bandwidth  sometimes  calculated  during  ducting  conditions  should  be 
viewed  with  caution  before  applying  it  to  an  error  rate  prediction  curve 
because  the  fade  structure  is  usually  very  shallow  and  seldom  plunges 
into  the  noise.  The  confuted  cross  correlation  coefficients  so  obtained 
are  normalized  by  subtracting  out  the  strong  dc  term  (average  signal 
strength)  in  computation.  Thus  a  relatively  meaningless  cross  corre¬ 
lation  coefficient  might  be  calculated  that  cannot  be  applied  to  the 
prediction  of  errors  when  such  error  predictions  are  based  on  Rayleigh 
fading  which  is  decidedly  not  the  case  during  ducting. 
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APPENDIX  A 

THE  BELLO  MODEL  FOR  THE  TROPOSCATTER  CHANNEL 
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APPENDIX  A 


THE  BELLO  MODEL  FOR  THE  TROPOSCATTER  CHANNEL 

A  fairly  recent  addition  to  the  literature  concerning  troposcatter 
channels  has  been  made  by  Bello.  As  this  model  was  not  treated  in  detail 
in  the  prior  interim  reports  of  this  contract,  it  will  be  considered  here, 
especially  as  it  pertains  to  the  four  links  for  which  experimental  data 
have  been  taken. 

This  model  rests  on  the  useful  concept  of  the  delay  power  spectrum 
Q ( £)  of  a  troposcatter  channel  which  is  the  time  domain  response  of  a 
(linear)  channel  to  a  unit  impulse.  It  is  possible  to  derive  from  first 
principles  the  delay  power  spectrum  expected;  essentially  the  Booker-Gordon 
model  is  used  to  find  the  scattered  power  from  each  shell  of  constant  time 
delay  (£).  These  surfaces  are  prolate  spheroids  with  transmitter  and 
receiver  at  the  foci.  The  delay  power  spectrum  is  analogous  to  the  notion 
of  multipath  spread  and  in  fact  t]he  rms  value  of  Q  (£)  could  be  considered 
as  just  A  .  Most  importantly,  Bello  shows  that  under  typical  conditions 
expected  to  be  present  for  the  troposcatter  medium,  the  frequency  correla¬ 
tion  function  is  precisely  the  Fourier  transform  of  the  delay  power  spec¬ 
trum.  Thus  emphasis  is  placed  on  obtaining  Q  (£ )  since  it  is  a  simple 
extension  to  obtain  R(fi) ,  the  correlation  function. 

The  power  received  via  a  troposcatter  link  from  a  single  scatterer  will 
be  given  by  the  bistatic  radar  equation,  i.e., 

dp  -  2JL  «  <»-H 

R  T 

where  R  and  T  are  the  distances  from  the  transmitter  and  receiver  respective¬ 
ly  to  a  Booker-Gordon  "blob"  having  scattering  cross  section  a.  G  and  H  are 
the  two  antenna  gains. 

The  notation  used  in  this  model  is  indicated  in  Figure  A-l.  The 
scattering  cross  section  cr  is  a  sensitive  function  of  the  scatter  angle,  0, 
and  also  depends  on  the  air  density  in  the  common  volume  or  the  height  h 
above  ground,  and  the  volume  dV  of  the  "blob."  Bello  writes 

<7  (9)  =  —  (A-2) 

0^ 

where  the  exponent,  m,  obtained  experimentally,  has  been  assumed  to  be 
approximately  five.  The  inverse  h  dependence  is  also  empirical,  although 
Bello  states  that  the  exact  form  is  fairly  insensitive. 
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Simplifying  to  a  cylindrical  co-ordinate  system  (i.e.,  cylindrical 
ellipsoids  in  lieu  of  prolate  spheroids)  ,  a  differential  volume  is  found 
to  be 

dV  ^  - — - - —  d^  df  dz  (A-3) 

2  i/4>y 

where  8  is  a  simplifying  notation  corresponding  to  a  normalized  path  delay: 


6  =  (c  =  vel.  of  light) 

D/C 

~  t± 

~  2  ‘ 

Using  this  notation,  R  and  T  are  approximately 


R  ”  °  !>  +  2  t/t 


D 


-ii/i _ 

*+  2j/V 


(A-4) 


(A-5) 

(A-6) 


Finally,  the  scattering  angle  0  is  very  nearly 

0  2r  ^  +  26/i  (A— 7) 

Thus,  the  power  scattered  toward  the  receiver  from  a  single  blob  of 
volume  dV  is 


dP 


G(Mto)  H(<M>o) 

s2  em~2 


d'i'  d{  dz 

T 


(A-8) 


The  quantities  ^o  and  ^o  do  not  correspond  exactly  to  the  normal  usage  of 
the  term  takeoff  angle.  As  seen  in  Figure  A-l,  $  and  <t>  are  measured 
from  the  straight  line  D  rather  than  up  or  down  from  the  smooth  earth 
tangent  at  the  two  foci. 


The  total  power  reaching  the  receiver  from  a  shell  characterized  by  £ 
is  then  equal  to  the  value  of  Q<£)  for  £  to  £+d£  and  is  obtained  by  inte¬ 
gration  over  z  for  some  reasonable  distance  and  over  ^ .  That  is, 

h 


Q({)  d« 


G(yfr-^o)  HC^-fo) 


m-2  , 

e  *l> 


dt fr  d{ 


(A-9) 


o  and  ^1  correspond  to  the  two  extremes  of  for  possible  scatter  toward 
the  receiver.  Below  't'o  no  power  can  propagate  past  the  transmitter  horizon; 
angles  above  ^1  correspond  to  a  receiver  angle  less  than  ^ o .  In  fact,  from 
the  above  A-4. 
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(A-10) 


*1 


Expressing  and  0  in  terms  of 


(from  A-4,7)  Q  ({)  becomes 


Q({)  =  $(<5)  ~ 


G (0-^o)  H(2S/*  -to) 


dtji 


(A-ll) 


For  calculational  convenience,  the  dummy  ^  may  be  replaced  by  x  =  25 

yielding 


§(8)  = 


g (l+m/2) 


f \fi. 


x/TJ/^c 


G(xyTT  -  ^o)  h cy/iT/x  -<M 


o 

V^ 


(x  +  -) 

X 


TTF2" 


dx  (A-12) 


Bello  assumes  a  symmetric  Gaussian  shape  for  the  two  beam  patterns  so  that 


G(f)  =  H(<|»)  *  exp (-4>2/o2) 


and 


a  ~  0. 6 (HPBW) 

The  above  expression  is  not  simply  integrable;  however,  numerical  inte¬ 
gration  is  straightforward.  Specifically,  for  the  results  that  follow, 

(XTRAN)  four  point  Gaussian  quadrature  was  used  with  Q(6)  found  for  512 
evenly  spaced  samples  of  8  .  The  range  of  6  was  chosen  from  (~'/t><fto/2)  to 
a  point  where  Q(8)  was  down  by  a  factor  of  109  from  its  peak  value.  Care 
must  be  exercised  to  insure  that  the  integral ' s  upper  limit  is  always 
greater  than  or  equal  to  the  lower  limit.  Q(5)  thus  found  will  not  be 
normalized  but  in  general  only  its  shape  is  important. 

Of  prime  importance  is  not  Q ( 5)  but  R(ft)  which  is  the  Fourier  transform 
of  Q.  This  is  again  a  simple  process  (using  a  computer)  with  the  Cooley- 
Tukey  Fast  Fourier  Transform  algorithm  employed  here.  The  entire  integration- 
transform  program  is  not  listed  in  this  report  but  is  available. 

For  test  links  of  this  program,  the  predicted  correlation  functions 
using  the  Bello  model  are  shown  in  Figure  A-2.  The  takeoff  angles  were 
calculated  assuming  4/3  earth  radius  and  not  considered  identically  zero. 

It  is  obvious  that  the  predicted  bandwidths  are  considerably  larger  than  the 
median  correlation  bandwidths  measured  here.  To  more  accurately  fit  the  data, 
several  modifications  might  be  attempted.  The  beamwidths  used  were  the  actual 
beamwidths  and  these  could  be  reduced;  however,  if  the  effective  beamwidths 
are  used,  the  multipath  spread  will  decrease  and  this  will  broaden  rather 
than  decrease  the  correlation  bandwidth.  The  parameter  m  might  also  be  altered, 
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especially  since  it  is  empirical.  Figure  A-3  displays  the  effect  of  various 
values  of  m  for  one  link,  Point  Petre,  C-band.  For  m  equal  to  3  or  1,  the 
situation  is  improved  although  still  wider  than  the  median  bandwidths 
observed. 

It  should  be  noted  that  these  findings  are  in  agreement  with  another 
comparison  of  the  Bello  model  with  a  previously  conducted  test  link, 

Ontario  Center-Youngstown,  (85  mile)  L-band  using  28  foot  antennas^,  in 
that  series  of  tests,  the  predicted  correlation  bandwidth  was  also  found 
considerably  larger  than  that  measured.  Since  the  Bello  model  gave  an 
excellent  fit  to  a  longer  L-band  path  (170  mile)  in  that  same  series  of 
tests,  it  probably  can  be  made  to  fit  the  C-and-X  band  data  in  this  re¬ 
port  with  some  modification.  One  relation  that  was  not  investigated  was 
the  1/h  dependence  of  <r.  This  dependence  might  not  be  as  insensitive 
as  thought. 
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Figure  A-2.  Frequency  Correlation  Function  Using  the  Bello  Model  for  the 
Four  Test  Links,  X-  and  C-Band 


1.0 


.8 


.6 


.4 


.2 


Figure  A-3. 


MHz 

Frequency  Separation 

Frequency  Correlation  Function  Using  the  Bello  Model  for 
Point  Petre,  C-Band  and  Varying  the  Parameter  m 
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cluud  ceiling  of  unknown 
height 


WLATHtRCOH  MS 


*  1  or  nidt  > 

f  Ihundcruorm 
()  Squill 
K  Riin 
R  W  Riin  showers 
7 R  l  ree/ing  rjm 
1.  Drizzle 
/L  I  rec/mg  drizzle 
S  Snow 
SP  Snow  pcflets 
1C  Ice  crystals 
SW  Snow  showers 
Si  i  Snow  grains 
P  Steel 
A  Hill 
AP  SfTlill  hill 

P  Pog 
IP  Ice  log 

CiP  C<  round  fog 
BI)  Blowing  duM 
BN  Blowing  sand 
BS  Blowing  snow 
BY  Blowing  spray 
K  Smoke 
H  Hi/e 
D  Dust 


WIND  COLUMNS 

Directions  arc  those  from 
which  the  wind  blows,  indi¬ 
cated  m  tens  of  degrees 
from  true  North;  i,  e..  09 
for  Ea«t,  IB  for  South,  27 
for  West  Entry  of  00  in 
the  direction  column  indi¬ 
cates  calm. 

Speed  is  expressed  in  knots, 
multiply  by  I  IS  to  convert 
to  miles  per  hour 


ADDITIONAL  DATA 

Other  observational  data  contained  in  records  on  file  can  be  furnished  at  cost  via  microfilm  or  microfiche 
copies  of  the  original  records  Inquiries  as  to  availability  and  costs  should  he  addressed  to 
Director.  National  Weather  Records  Center,  Pederal  Building.  Asheville.  N  C  28801 
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LOCAL  CLIMATOLOGICAL  DATA 

U.S.  DEPARTMENT  OF  COMMERCE 

MAURICE  H.  STANS,  S«cr«tory 

ENVIRONMENTAL  SCIENCE  SERVICES  ADMINISTRATION 

environmental  data  service 
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*  Extreme  tem|«-rnturi‘s  for  the  month  Mav  In*  the  la*l 
of  more  than  one  occurrence 
-  Below  zero  temperature  or  negative  depauure  front 
normal. 

J  ,  70  at  Alaskan  stations. 

+  Also  on  an  earlier  date,  or  date- 
X  Heavy  fog  restricts  visibility  t-  >4  mile  or  !.•*»**. 

T  In  the  Hourly  Precipitation  table  and  in  column* 
9,  10.  and  11  indicates  an  amount  too  small  to 
measure. 

The  season  for  degree  days  begins  with  Julv  for  heating 
and  with  January  for  cooling. 

Data  in  column*  6,  12.  13,  14,  and  lf»  are  ba-ed  on  8 
observations  per  day  at  3-hour  intervals. 

Wind  directions  are  those  from  which  the  wind  blows. 
Resultant  wind  Is  the  vector  sum  of  wind  directions 
and  speeds  divided  by  the  number  of  observations. 
Figures  for  directions  are  tens  of  degrees  from  true 
North  ;  i.e.,  09  ~  East,  18  -  South,  27  •  West,  30  North, 
and  00  =  Calm.  When  directions  are  in  tens  of  degri-es 
in  Col.  17,  entries  in  Oil.  K>  are  fastest  observed 
1-minute  speeds.  If  the  /  appears  in  Oil.  17,  speeds 
are  gusts. 
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Am  error-  detected  will  lie  corrected  and  change.-  ;n 
summary  data  will  he  annotated  in  the  annual  -umnai-y. 


^Subscription  Price  Lorn  I  Climatological  Daln  S 1  o  • 

|  per  vear  including  anminl  Summary  if  published  j 
;  Singlp  copy:  in  cents  for  monthly  Summary;  If* 
j  cents  for  annual  Summary.  Check-  or  mon-v  orders 
should  Ik-  made  payable  and  remittances  and  n  rre  j 
spojidence  should  Ik*  sent  to  the  Superintenilenl  of 
|  Documents.  1'  S  finvermnent  Prin*  mg  OflW 
1  Washington,  1>  J  20HVJ. _ J 


]  certify  that  this  is  an  official  publication  »f  the 
Environmental  Science  Services  Administration,  and  is 
compiled  fmm  records  on  fib*  at  the  National  Weather 
Records  Center.  Asheville,  North  Carolina  28801. 

•/J. 

Director,  National  Weather  Record-  Center 
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•  Extreme  temper*  lures  for  the  month  Mav  tie  the  last 
of  more  than  on*  occurrence 

-  Below  zero  temperature  or  negative  departure  from 
normal. 

2  >  70°  at  Alaskan  stations. 

+  Alao  on  an  earlier  date,  or  dati-a. 

X  Heavy  fog  restricts  visibility  to  <4  mile  or  less 
T  In  the  Hourly  Precipitation  table  and  in  columns 
9,  10,  and  11  indicates  an  amount  too  small  to 
measure. 

The  season  for  degree  days  begins  with  July  for  heating 
and  with  January  for  cooling. 

Data  in  coiumna  8,  12,  IS,  14.  and  15  are  based  on  8 
obaervationa  per  day  at  3-hour  intervals. 

Wind  direction*  are  those  from  which  the  wind  blows. 
Resultant  wind  Is  the  vector  sum  of  wind  directions 
and  speeds  divided  by  the  number  of  observations. 
Ftfurea  for  directions  an*  tens  of  degrees  from  true 
North ;  i.e.,09  =  East,  18  =  South,  27  -=  West.  36  -  North, 
and  00  =  Calm,  When  directions  are  in  tens  of  degrees 
in  Col.  17,  entries  in  Col.  16  are  fastest  observed 
1-minute  speeds.  If  the  /  appears  in  Col.  17.  speeds 
are  gusts. 


Any  errors  detected  will  be  corrected  and  changes  in 
summary  data  will  be  annotated  in  the  annual  summary. 


SUMMARY  BY  HOURS 


Subscription  Price:  Local  Climatological  Data  81.00] 
per  year  including  annual  Summary  if  published. 
Single  copy:  10  cents  for  monthly  Summary;  15j 
cents  for  annual  Summary.  Checks  or  money  orders) 
should  be  made  payable  and  remittances  and  corre 
apondence  should  be  sent  to  the  Superintendent  of) 
Documents.  U.  S.  Government  Printing  Office, 
Washington.  D.  C.  20402. _ _ _ 


I  certify  that  this  is  ««  official  publication  of  the 
Environmental  Science  Services  Administration,  and  is 
compiled  from  records  on  file  at  the  National  Weather 
Records  Center,  Asheville,  North  Carolina  28801. 

Iv'Mu—  / /. 

Director,  National  Weather  Records  Center 
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Washington,  I>.  C.  20402. 


I  certify  that  this  is  an  official  putilir.it i»»n  of  the 
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compiled  from  records  on  file  at  the  National  Weather 
Records  Center.  Asheville.  North  Carolina  28801 

kfUAi —  ;*/ 
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Also  on  an  earlier  date,  or  date*. 

X  Heavy  fop  restricts  visibility  to  14  mile  <*r  less 
T  In  the  Hourly  Precipitation  table  and  in  columns 
9.  10.  and  11  indicates  an  amount  -mall  to 
measure. 

The  season  for  degree  days  begins  with  Julv  for  heating 
and  with  January  for  cooling. 

Data  in  columns  fi.  1  13,  14.  and  15  are  ha-ed  on  S 

observations  per  day  at  3-hour  intervals 
Wind  directions  are  those  from  which  the  wind  Mows, 
Resultant  wind  is  the  vector  sum  of  wind  direction- 
and  speeds  divided  by  the  number  of  observations 
Figure*  for  directions  are  lens  of  degree-  from  true 
North ;  i.e..  09  =  East,  18  =  South,  27  West.  3fi  North, 
and  00  -  Calm.  When  directions  are  in  tens  of  degrees 
in  Col,  17,  entries  in  Col.  16  an*  fastest  observed 
l*minute  speeds.  If  the  /  appears  in  Col.  17.  speeds 
are  fuats. 


Any  error-  detected  will  be  corrected  and  changes  in 
summary  data  will  l»e  annotated  in  the  annual  summary 


[Subscription  Price:  I^k-hI  Climatological  Data  $1.00 
per  veat  including  annual  Summary  if  published. 
Single  copy  10  icnts  for  monthly  Summary  15 
cent-*  for  annual  Summary.  Checks  or  money  orders 
should  1m>  made  payable  and  remittances  and  c«»rre 
spondenrr  should  tie  sent  to  the  Superintendent  of 
Documents,  f  S.  Government  Printing  office. 
Washington.  1)  C.  20402.  _ _  .  __  . 


I  certify  that  this  is  an  official  publication  of  the 
Environmental  Science  Services  Administration,  and  is 
compiled  from  records  on  file  at  the  National  Weather 
Records  Center.  Asheville.  North  Carolina  2W*ni 

f/. 

Director.  National  Weather  Records  Center 
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ADDITIONAL  OAT  A 

Other  observational  data  contained  in  record'  on  file  can  he  turntshed  at  emt  via  microfilm  or  rmcroliche 
copies  of  the  original  records  Inquiries  a*  to  avdilanilnv  and  costs  should  he  addressed  u* 

Director.  National  VX  cat  her  Records  l  enter,  federal  Building.  Asheville.  N  ('  2HH0 1 

STATION:  ROCHESTER  N  Y  YEAR  k  MONTH:  70  01 
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APPENDIX  C 

EXTENDED  CORRELATION  BANDWIDTH  PROGRAM 
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MAIN  PROGRAM:  TROPO 


Inputs : 

ALTW= altitude  in  meters  above  sea  level  of  western  most  station 
ALTE=altitude  in  meters  above  sea  level  of  eastern  most  station 
TF= temperature  in  degrees  F. 

H=hjmidity  in  percent 

P=atmospheric  pressure  in  inches  mercury 
BW=antenna  beam  width  in  degrees 

SF=scale  factor  related  to  condition  of  the  troposphere 

Calculated  in  subroutines  or  input: 

S=distance  in  kilometers  between  sites  on  great  circle  path 
GAMMAW=  angle  in  radians  between  tangent  ray  and  lower  ray  over 
obstacles  for  western  most  site 
GAMMAE= angle  in  radians  between  tangent  ray  and  lower  ray  over 
obstacles  for  eastern  most  site 


Outputs : 

NS=surface  refractivity 

ALPHA=effective  antenna  beam  width  in  radians 
SUPPR=upper  ray  length  in  kilometers 
5LOWR=lower  ray  length  in  kilometers 
AC=multipath  spread  in  microseconds 
THETA=scatter  angle  in  radians 
RE=envelope  correlation  function 
F=frequency  separation  in  MHz 

SUBROUTINE:  ARC 


Inputs : 

Longitude  and  latitude  in  degrees,  minutes,  and  seconds  of  western  most  site 
Longitude  and  latitude  in  degrees,  minutes,  and  seconds  of  eastern  most  site 
Zl ,  Z2=codes  for  sites  east  or  west  of  Greenwich 

M=increment  in  miles  along  great  circle  path  for  print  out  of  bearing  and 
distance 

Outputs : 

S=distance  in  kilometers  between  sites  on  great  circle  path 
B7=bearing  east  in  degrees  for  great  circle  path 
C3=bearing  west  in  degrees  for  great  circle  path 

SUBROUTINE :  ANGLE 


Inputs : 

OBHT=height  of  obstacle  in  meters  above  sea  level 
OBDIST=distance  in  meters  from  site  to  obstacle 

Outputs : 

GAMMAW=angle  in  radians  between  tangent  ray  and  lower  ray  over  obstacles 
for  western  most  site 

GAMMAE=angle  in  radians  between  tangent  ray  and  lower  ray  over  obstacles 
for  eastern  most  site 
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